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Abstract 

We show that heavy pure states of gravity can appear to be mixed states to 
almost ah probes. For AdSs Schwarzschild black holes, our arguments are made 
using the field theory dual to string theory in such spacetimes. Our results 
follow from applying information theoretic notions to field theory operators 
capable of describing very heavy states in gravity. For half-BPS states of the 
theory which are incipient black holes, our account is exact: typical microstates 
are described in gravity by a spacetime "foam", the precise details of which 
are almost invisible to almost all probes. We show that universal low-energy 
effective description of a foam of given global charges is via certain singular 
spacetime geometries. When one of the specified charges is the number of D- 
branes, the effective singular geometry is the half-BPS "superstar" . We propose 
this as the general mechanism by which the effective thermodynamic character 
of gravity emerges. 
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There are twenty-five orthographic symbols. That discovery enabled mankind, 
three hundred years ago, to formulate a general theory of the Library and thereby 
satisfactorily solve the riddle that no conjecture had been able to divine — the 
formless and chaotic nature of virtually all books. . . . The Library is total and 
its shelves register all the possible combinations of the twenty-odd orthographic 
symbols (a number which, though extremely vast, is not infinite). 

— Jorge Luis Borges, in "The Library of Babel" 

1 Introduction 

The microcanonical accounting of the entropy of black holes posits an enormous 
number of degenerate microscopic states. In terms of these states, the information loss 
paradox may be evaded by showing that a pure initial state collapses to a particular 
pure black hole microstate whose exact structure can be deduced by suitably subtle 
measurements. Given the importance of the problem, it is crucial to ask: What 
do pure microstates look like and what sorts of measurements can distinguish them 
from each other? Here we report that in string theory almost no probes are able 
to differentiate the microstates of a black hole. Thus, if spacetime is placed in a 
very heavy pure state, it will appear mixed — i.e. like a black hole — to almost all 
finite precision measurements. This explains how the existence of pure underlying 
microstates and the absence of fundamental information loss are consistent with the 
semiclassical observation of a thermodynamic character of black holes. 

In Section 2, we analyze so-called large Schwarzschild black holes in AdS space- 
times. These black holes have horizon size bigger than the scale set by the AdS 
curvature. They are known to come into equilibrium with their thermal radiation 
because AdS geometries create an effective confining potential pp. Thus such black 
holes are stable, and we can ask how the underlying states that give rise to the large 
entropy can be identified via quantum mechanical probes. We will examine black 
holes in five dimensions because string theory on AdSs x admits a dual description 
in terms of a superconformal SU{N) gauge field theory propagating on a three-sphere 
12] . Here N is related to the AdS curvature scale 1/i via the string coupling gs and the 
string length ig as i'^ = 27rgsNig. When the length scale i is large, which we require, 
so too is A^. In the dual field theory, we characterize the structure of operators of 
conformal dimension 0{N'^) that create black hole microstates. Almost all such op- 
erators belong to a "typical set" , which is characterized by the statistical randomness 
of the polynomial in elementary fields that defines the operator. Sanov's theorem 
from information theory shows that deviations from typicality are exponentially sup- 
pressed in the conformal dimension. Atypical operators are vanishingly improbable 
in the large- limit. We then argue that almost no probe correlation functions can 
distinguish two typical states from each other. We also explore what questions about 
black holes are easy to answer using probe measurements. 
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Our arguments about AdS-Schwarzschild black holes are made in the dual field 
theory because the underlying non-supersymmetric microstates have not been given 
a geometric description. For this reason, in Section 3 we examine half-BPS states 
of the theory which have been explicitly described in both the AdS spacetime and 
the dual field theory. Very heavy half-BPS states are incipient black holes — a small 
amount of energy beyond the BPS limit is expected to create a horizon. In field 
theory, these states are described in terms of free fermions in a harmonic potential, 
or equivalently, in terms of Young diagrams enumerating gauge-invariant operators 
[31 Hj. Like the Schwarzschild case, almost all states can be characterized by the 
statistical randomness of the operators that make them. Specifically, almost all such 
states are associated to Young diagrams whose shapes are small fiuctuations around a 
certain limiting curve. This characterization of typical states is done in two different 
statistical ensembles: with and without a restriction on the number of columns in the 
Young diagrams. Physically, the second ensemble corresponds to fixing the number of 
D-branes in the state. The infinite temperature limit in the second ensemble describes 
the set of typical states underlying the superstar [S]. 

In Section 4 we relate the exact quantum mechanical description of the mi- 
crostates, and their averages studied in Section 3 with their semiclassical gravitational 
description.^ We develop techniques for characterizing properties of typical states in 
terms of phase space distributions and then establish a dictionary with the classical 
asymptotically AdS geometries to which they correspond. We establish the criteria for 
half-BPS states to have such semiclassical descriptions and show that a typical very 
heavy state corresponds to a spacetime "foam". The precise structure of the foam, 
and whether it looks more quantum or classical, depends on the observable. We argue 
that almost no semiclassical probes will be able to distinguish different foam states, 
and that the resulting effective description gives a singular geometry. For spacetimes 
in which the overall mass and charge are specified the effective geometry is a new 
singular configuration that we call the "hyperstar" (see also [E]). For spacetimes in 
which the number of D-branes is additionally specified, the effective geometry is again 
singular, giving rise in the maximum entropy limit to the well-known "superstar" j3] 
geometry. 

In Section 5, we argue for the universality of a subset of observables in the one- 
half BPS sector. We also discuss particular observables which are more efficient 
at discriminating between different microstates, but argue against the feasibility of 

^Our work is related to the proposals of Mathur and others [HJ |7| HJ |H1 El ^2 El in the context 
of the microstates of the D1-D5 system. These works demonstrated the existence of non-singular 
geometries without horizons matching the quantum numbers of certain D1-D5 microstates. The work 
of El '^^^ be regarded as providing the analogous half-BPS geometries in AdSs x S^. Our perspective 
on these results is that they provide a classical configuration space that requires quantization. We 
are using the field theory dual to AdSs x to carry out this quantization and carefully study the 
semiclassical limit. 
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constructing such probes. 

This universal low-energy effective description of underlying smooth quantum 
states is the origin of the thermodynamic character of heavy gravitating systems. 



2 Schwarzschild black holes 



Consider the metric for a Schwarzschild black hole in AdSt^: 



1 + ^-^ 



1 + — - 

£2 ^2 



dr^ + r^dn^ + d{S^i)), 



where i is the radius of AdS^; and of S^. The horizon is at 



V2 



1 + 11 + ^ 



2\ 1/2' 



1 1/2 



The difference between mass of the black hole and the mass of empty AdSs is 

M 



(2) 



(3) 



Dropping numerical constants, the five-dimensional Newton constant is Gs ~ g\c^'' . 
Also recall the AdS / CFT dictionary relating string theory on AdSs to A/" = 4, SU (N) 
super- Yang-Mills theory at large- A^: 

g^N ~ Qym^ = a = fixed 't Hooft coupling, 

i ^ U9sNy/\ (4) 
Mi ~ A = conformal dimension of the dual operator. 

For a large black hole, c:^ i which implies tq ~ i. (For a small black hole, tq ~ 
Th -C i] such small black holes are unstable to localization on the sphere and we will 
disregard them here.) If a spacetime of this mass is described by a pure state in the 
CFT, the conformal dimension of the operator making the state is 



A 



M£ ^ — 



9i 



(5) 



which is independent of the Yang-Mills coupling. In view of the fact that these 
states are not supersymmetric, and should renormalize, it is remarkable that (jH)) is 
independent of the Yang-Mills coupling. 

From the CFT point of view, a pure underlying black hole microstate is described 

as 



Imicrostate) = C»|0); A(C>) ~ A^l 



(6) 



It is worth comparing this to the conformal dimensions of various more familiar states 
in the AdS/CFT correspondence: 
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Supergravity states in AdS^+i 




Small strings 


E ~ 1/4 — ^ A ~ ^/^s ^ {gsNf^ ~ A^/^ 


Small black holes (unstable) 


M ~ t'l — ^ A ~ '''^ ~ Ar2;o_ 


D-brane states (giant gravitons) 


A ~ A^. 



Large black holes are extremely heavy states, with a conformal dimension that grows 
in the large- limit like N D-branes. The black hole entropy is 



5 ~ ^ ~ ^ ~ ^^Gio ~ iV' ~ A, (7) 

SO we expect that the field theory dual to AdSs space has 0{e^^) operators of confor- 
mal dimension A = N"^ that are neutral under the global symmetries of the Yang- Mills 
theory. 

2.1 Typical states of Schwarzschild black holes 

The fields of = 4 SYM are the gauge fields A^j,, three complex adjoint scalars 
X, Y, Z, and four complex Weyl fermions ipa- Gauge invariant operators can be con- 
structed as products of traces of polynomials built out of the fields of the theory and 
their derivatives, with Lorentz and gauge indices appropriately contracted. We wish 
to construct operators of dimension A^^ that are not charged under either the 5*0(4) 
rotation group or the 5*0(6) R-symmetry. An example gauge invariant operator con- 
structed from the scalars is 

o = (tr[A:A:^yy^z])(tr[A:^y^y^yz^A:xy])(tr[xA:^A:^])(tr[y^yy^zyz^A:^]) .... 

(8) 

Derivatives could be included also as well as fermions and gauge fields, but we ignore 
these for the moment, since our arguments will not depend on the details of particular 

operators. 

The operators in question are not BPS since the Schwarzschild black hole breaks 
supersymmetry; thus terms like XX^ may appear. (Since R-charges effectively count 
differences in say the number of Xs and Xh, a chargeless operator should have equal 
numbers of a field and its conjugate.) In the absence of supersymmetry, operator 
dimensions renormalize, and thus the classical dimension could be different from the 
full quantum corrected dimension A = N'^ that we are seeking. However, it is possible 
that for such heavy operators renormalization is minimal. One can argue for this in 
three steps: 

1. The conformal dimensional A = A"^ is independent of the couphng. Since this 
is derived from the mass of the dual black hole in the semi classical limit, we 
learn that at least at large 't Hooft coupling the conformal dimension does not 
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receive corrections as the coupling is changed. This leaves open the possibility 
that the renormalization of such operators is large at weak coupling. 

2. Large black holes in AdS, such as the ones we study, are usually modeled as 
thermal ensembles in the dual field theory. It is known that the thermal free 
energy only renormalizes by a factor of 3/4 from weak to strong ^H]- This 
suggests that the energies of the heavy states entering the thermal ensemble 
also do not renormalize substantially from weak to strong coupling. 

3. A possible explanation for this lack of renormalization despite the absence of 
supersymmetry is that such enormously heavy states are effectively classical 
because of cancellations between the contributions to the self-energy from ef- 
fectively random interactions between various pieces of large operators like (jH)). 

In view of these arguments we will start in free field theory and ask how to construct 
neutral operators of large conformal dimension. In any case, we will see that the basic 
result, that almost no probes can identify a state created by a very large operator, 
will not depend very much on the dimension being precisely A^^, and that interactions 
will not modify our conclusions. 

At large- A^, low dimension operators that are polynomials in traces are approxi- 
mately orthogonal in the sense that the states they create have small overlap. How- 
ever, at large conformal dimension the trace basis mixes heavily even in the free 
field limit. For example, an orthogonal basis of "giant graviton" D-brane operators 
of dimension ~ A^ is constructed from determinants and sub-determinants of the 
scalar fields [121 CZl dl UHj; even though these can of course be expanded as linear 
combinations of traces. In addition, in an SU{N) theory a trace of a product of 
more than A^ fields decomposes into a sum of smaller traces. We will not discuss the 
problem of finding an orthogonal basis of non-supersymmetric operators of dimension 
0(A^^) because this will not be necessary. We will argue that almost all operators 
of high dimensions have the properties we seek, and thus, in free field theory, any 
linear combination of them will do. Subject to these caveats, which we will return to 
later, each operator can be thought of as a sentence built of words (traces), each of 
which is a string of letters in the alphabet provided by the elementary fields of the 
theory and their derivatives. Black hole microstates are created by operators with 
very large dimensions (A ~ A^^) and hence (ignoring traces of gauge indices) involve 
polynomials of length of order A^^. Theorems in information theory characterize the 
structure of such polynomials. As the conformal dimension becomes large, almost all 
operators will belong to a "typical set" characterized by the statistical randomness 
of the sequence of fields in the polynomial. Specifically, for operators in the typical 
set the probability that a randomly chosen letter in the polynomial is x is given by 
the uniform distribution q{x) = 1/k, where k is the size of the alphabet. (For the 
moment we are ignoring derivatives as "letters in the alphabet".) In particular, this 
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ensures that for operators in the typical set, the number of Xs will be the same as 
the number of X'^s, etc. so the R-charges will automatically vanish. Deviations from 
typicality are exponentially suppressed. For example, the likelihood that a randomly 
chosen operator of length A has a letter distribution p{x) is given by 

P{p{x)) = 2-^^(Pll5) (9) 

where D{p\\q) is the Kullback-Liebler distance^ or relative entropy between probabil- 
ity distributions 

^(pik) = Ep(^)i°g2(^)- (10) 

This is an approximation at large A of the exact answer 

= ( Ap(l) Ap(2) (11) 

In the present case, the precise result requires some modification to include the 
distribution of traces and derivatives, and the fact that the different Yang- Mills fields 
have different conformal weights. For example, higher derivatives are permitted, thus 
enlarging the "alphabet", although equations of motion will relate many of these 
operators. A proper accounting of these effects is necessary in order to obtain e^^ 
polynomials of length N"^ that are necessary to give the black hole entropy. This 
counting was studied in detail in [2011^1221^ and by Legendre transforming the result 
for the free energy at high temperature and zero coupling we obtain the asymptotic 
density of states 

Ar(A) ^ et"^*^\ (12) 

For conformal weights A ~ this indeed yields the right e^^ degeneracy of opera- 
tors. 

Roughly speaking, the typical black hole microstate operator will look like a ran- 
dom polynomial with traces randomly distributed in it. Operators such as 

O' = {tr[XYXYXYXYXYXY . . .]){ti[XYXYXYXYXYXY ...])... 

with discernible order are exponentially rare at large conformal dimension. To gain 
some intuition for these statements consider single trace operators of dimension n 
built out of k species of scalar fields each having conformal dimension 1. Then the 
distribution of letters in the typical state is q{i) = 1/k. A randomly chosen operator 
will have a particular letter distribution p{i) = 1/k + 6i/n. (The 6i satisfy the con- 
straint Yli '^4 = since the word length is fixed.) Then Sanov's theorem j2S| measures 

^Calling this a distance is a slight abuse of terminology. is not symmetric and does not 

satisfy the triangle inequality, and so is not a metric distance. 

•^In |2J abelian theories are studied, whereas HOI 121] consider the nonabelian situation. 
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the probability of having a large deviation from the expected distribution q{x):'^ 



P{p{x)) = 2-"^(fll'') ~ 




Large deviations are therefore exponentially suppressed. We are simply describing a 
generalization of the central limit theorem — there are 0(^^~rifk) variations in the 
frequencies of different fields appearing in the long operator. This may appear large, 
but of course the ratio of the standard deviation to the mean, 1 / ^Jn/k, decreases as 
n increases. 

2.2 Probing a typical state 

Consider states \0) = O\0) created by operators of very high conformal dimension 
and a collection of probe operators Op. 

Claim: Correlation functions of probe operators Op computed in very heavy states 
\0) {i.e. {0\0p • • ■ Op\0)) depend only on the conformal dimension A and the global 
charges of Op and O, up to corrections of 0(e~^), for almost all states \0) and probes 
Op. 

In Section 5, we will explore this claim for half-BPS states after showing in Sections 
3 and 4 that the resulting effective description of underlying smooth quantum states 
in low-energy gravity can be a singular spacetime. Here, in view of the challenges of 
rigorously enumerating the microstates of a Schwarzschild black hole we will content 
ourselves with an argument rather than a proof. 

For purposes of argument, consider the two point correlation function of a light 
probe Op = tr[XX] in a heavy state O = tr[XYZXXYYZ ■ ■ ■]. In free field theory, 
Feynman diagrams contributing to the correlation function are computed via Wick 
contractions of fields in Op and O, e.g., 



(0| tr[XYZXXYYX • • f tr[XX]^ tr[XX] tr[XYZXXYYX ■ ■ ■] |0) (14) 

In free field theory each pattern of such Wick contractions makes a separate contri- 
bution to the overall correlation function. Contractions that are planar give contribu- 
tions that have the leading order in N and require the probe operator Op to appear 
as a sub-word within the state operator O as the minimum number of index lines are 
broken in this way. The contribution from planar diagrams is thus controlled by the 

^The reason why this probabiHty is not proportional to the length of the random polynomial is the 
normalization chosen for the fluctuations Si. If wc would have parametrized them as p(i) = l/k + 6i, 
this probability would indeed be manifestly suppressed by the length of the operator n, that is, by 
its conformal dimension A. 
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number of times the probe appears within the state, a quantity that is completely 
determined by the statistics of long random polynomials. In fact, at each order in 
the 1/N expansion the contribution of different patterns of Wick contractions is de- 
termined by how often each of these patterns occur. This is also controlled by the 
statistics of random polynomials. Adding perturbative interactions does not change 
the argument — one must then include the interaction vertices, but the number of 
times each Feynman diagram pattern occurs is still determined by statistics. In ef- 
fect the long polynomials making the state are creating a random "background" with 
which the probe interacts.^ The resulting statistical character of the computation 
guarantees that correlation functions in any typical state will take a universal value 
depending largely on the charges that control the overall frequencies of patterns in 
the state operator. Deviations from this universal value will be exponentially sup- 
pressed as they will depend on the tiny structural differences between specific heavy 
operators and the idealized typical state. 

The argument above was made in the context of single trace operators made from 
scalar fields, but it is clear that increasing the alphabet to include all the fields of 
the Yang-Mills theory, derivatives and randomly sprinkled traces will not change the 
essential statistical character of the computation. (Of course, the precise details of 
the Sanov theorem (jH)) will be modified since multiple derivatives will give rise to 
an alphabet size that grows without bound.) Neither will the computation of higher 
point correlators. The mixing that is expected between multi-trace operators of very 
high dimension only improves the argument for the following reason. We should 
expect all operators in the trace basis that have given quantum numbers to mix with 
about equal strength. Thus, in effect, the mixing will average correlation functions 
over the typical set of large dimension operators, further emphasizing the universality 
of the correlation functions. Interestingly, following ©, deviations from universality 
in the correlations should be controlled by the small parameter 



Since the entropy of the Schwarzschild black hole of mass A is S* ~ A^^, this sup- 
pression is exactly as expected for a large statistical system with e'^ underlying mi- 
crostates. 

Since our arguments for universality of correlation functions will work for any 
sufficiently heavy state, why are we describing a black hole and not just a thermal gas? 
In a conventional thermal gas, an individual probe with energies much larger than the 
thermal scale passes through the system unaffected. By contrast, heavier probes of 
black holes are increasingly likely to be absorbed. Indeed, a high conformal dimension 
probe does not decouple from a heavy state. Rather, the many fields entering into 

^This is reminiscent of computations of open string interactions in the presence of D-branes in 
Yang-Mills theory [n|[I3[Tn|. 
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the operator will individually interact with the fields in the state operator, increasing 
the effective strength of the interaction. This is simply telling us that the field theory 
knows about gravity, and that the effective coupling involves the energy of states. The 
standard field theory statement of decoupling, that high energy probes decouple from 
a low temperature thermal bath, applies to individual probes of very high momentum 
passing through a long- wavelength thermal background. This is a different notion of 
high energy than the one relevant here, namely the conformal dimension of a very 
heavy probe. 

To actually see the precise difference between a black hole and a thermal state 
such as a neutron star, and in particular the presence of a classical horizon, it will 
be necessary to investigate the behavior of the correlator as a function of Yang-Mills 
coupling. The reason for this is that any conventional star will collapse to become a 
black hole if the Newton constant is sufficiently large. Thus we should expect that as 
the string coupfing (which is dual to the Yang-Mills coupfing) is increased, standard 
thermal states can become black holes. While some sort of "phase transition" at finite 
coupling is presumably involved, it is possible that some of the classic properties of 
black holes, such as the difficulty of extracting information can already be understood 
as described above. This is because the universality of correlation functions is being 
governed by the randomness of the polynomials in fields that make up an operator 
along with the randomized distribution of traces that will produce an ensemble of 
small and large words in the complete operator. Indeed, at least in perturbation 
theory at finite coupling, it is clear that the statistical character of the arguments 
that we have made does not change. Perhaps all of this should be read as saying 
the difficulty of identifying the underlying state of a complex system is not by itself 
the signature of the presence of black hole with an "interior region" that is causally 
disconnected from infinity in the classical limit. 

2.3 Askable questions 

Given the statistical nature of the operators O creating black hole microstates, how 
much information can different kinds of probes give us about this state? Recall 
that the conformal dimension of (9 is A ~ N^. As argued above, the correlation 
functions of a probe can be used to piece apart the patterns of fields appearing 
within the definition of O. Supergravity probes with A ~ 0(1) will be completely 
ineffective at this task for several reasons. First, as we have argued, their correlators 
will be universal up to corrections of order e~^^ . Of course, one might compute the 
correlators of many different probes to separate the different microstates. However, 
simple counting shows that there are insufficient light operators of dimension 0{1) 
to distinguish the e^' expected microstates. Stringy probes (A ~ (gsNY^^) will not 
be much better since the number of possible patterns of letters within microstates of 
dimension ~ N'^ is so large. 
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Brane probes (A ~ A^) will be better placed to disentangle the microstate. In par- 
ticular the long stretches of fields appearing in the definition of the brane operator 
might or might not match sequences of fields in the microstate and these two condi- 
tions will lead to somewhat different correlation function responses. One might also 
use other black hole microstates {i.e. other operators of dimension 0(A^^)) as probes. 
The correlation functions of almost all such probes would involve the interaction of a 
random polynomial (the probe) with a different random polynomial (the state), lead- 
ing to a correlation function of universal form. However, if the experimenter happens 
to choose a probe that is identical to the black hole microstate, a huge response will 
result. Using black hole microstates as probes in this way is difficult because there 
should be e^^ such operators. Thus we see that it is far easier to detect what the 
state of the black hole is not rather than what it is.^ 

To learn what the state of a black hole is, we might make measurements with 
many different probes over a long period of time. This is in analogy with measuring 
higher moments to determine a function. We might also ask statistical questions like 
"How should a set of K probes be designed to optimally classify the set of black holes 
of mass M into K classes?". Our discussion has focused on typical microstates of 
fixed energy. It is easy to generalize to typical states in which other charges are also 
fixed. For example, fixing the R-charges to take non-zero values will constrain the 
number of X fields minus the number of X"^ fields appearing in the operator. For 
states with small R-charges, these can be regarded as introducing a small number of 
defects (a surplus of one of the fields) into a typical (therefore chargeless) random 
string in all possible ways. Of course, we could consider ensembles of operators in 
which more intricate patterns are fixed {e.g., one can consider operators in which 
the pattern XYYZ occurs twice as frequently as the pattern YXXZ.) It is possible 
that such restrictions, which are not related to specifying global charges, can be 
associated to the specification of quantities such as higher moments of fields in the 
dual gravitational description of black hole microstates as proposed by Mathur et al 

laiziiHii. 

Overall, the main lesson is that almost all probes yield almost no information 
about the structure of a black hole microstate. Thus, even if there is no fundamental 
information loss from an underlying pure heavy state in gravity, it will be extraor- 
dinarily difficult to identify what the state is. Indeed, the only information about 
the state operator that we can readily extract are conformal dimension, spin, and R- 
charge. On the gravity side, these correspond to mass, angular momenta, and global 

^If the underlying theory is actually integrable, there will be an infinite number of possibly 
nonlocal conserved quantities. Measuring these in detail would in principle identify the black hole 
state. However, even in this case, one would expect to have to make a very large number of very 
precise measurements to distinguish among microstates. From this perspective, our argument is 
simply that typical very heavy states are very similar to each other, and that very high precision 
will be needed to tell them apart. 
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charges. In an effective picture of probes with hmited energies making measurements 
in hmited regions of space and time, the standard classical picture of black hole with 
a causahy disconnected region should re-emerge. It is difficult to explicitly demon- 
strate all of this for Schwarzschild black holes because the underlying quantum states 
have not been constructed in spacetime and are difficult to analyze even in the dual 
field theory. Therefore, in Sections 3, 4 and 5 we will demonstrate how our consid- 
erations apply to half-BPS states of the theory which do not renormalize, and for 
which all computations, both in gravity and in the dual field theory, can be carried 
out explicitly. We will see that almost all heavy half-BPS states are described micro- 
scopically as a quantum "foam", whose universal effective description at low energies 
is a certain singular spacetime. The structure of the computations is identical to the 
discussion above which gives us some confidence that this is the correct picture even 
for Schwarzschild black holes. 



3 1/2-BPS States: Field theory 

Half-BPS multiplets of the A/" = 4 Yang-Mills theory transform in (0,p, 0) represen- 
tations of the SO (6) R-symmetry group. Highest weight states in each multiplet can 
be constructed as gauge- invariant polynomials in a complex scalar field X. Since the 
superfield X has conformal dimension A = 1 and carries R-charge J = 1 under a 
U{1) C ^0(6) each highest- weight half-BPS state has A = J. 

A convenient method for constructing an orthonormal basis of half-BPS operators 
utilizes an isomorphism between representations of U{N) and the symmetric group 
E]. It is an old observation due to Frobenius that the representations of the 
symmetric group are in one-to-one correspondence with the partitions of the integers. 
These are conveniently considered using Young (Ferrers) diagrams. Consider a Young 
diagram T with k boxes encoding one such representation of the symmetric group Sk- 
For example, for A; = 9 we might have 



1 


2 


3|4| 


5 


6 




7 


8 




9 







where we label the boxes of the Young diagram with integers 1, . . . , /c. We can con- 
struct an associated Yang-Mills operator 

OriX) = CrY: XB{a)Xl^XZ^^ ■ ■ ■ X^^ (17) 

a(T) 

where the sum runs over permutations of {1, ... , k} and Xr{(^) is the character (trace) 
of the permutation a in the representation R of the unitary group encoded by the 
same diagram T. Ot{X) is a Schur polynomial in the Xs. Such polynomials supply 
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a complete basis for the symmetric polynomials of degree kJ This is just the trace of 
X in the representation R. The coefficient Ct is chosen to normalize the two-point 
function so that (0|OT(a;)Or(2/)|0) = l/\x - y\'^^. All half-BPS operators can be 
written as (descendants of) linear combinations of the Ot- 

It can be shown that the constraint A = J implies that the half-BPS states can also 
be constructed as excitations of a Hermitian matrix in a harmonic oscillator potential 
13]. This can in turn be diagonalized to a theory of fermions {qi, . . . , g„} in a 
harmonic potential. The ground state of this system, the filled Fermi sea, consists of 
fermions with energies Ef = {i — l)huj + huj/2 for i = 1, . . . , N. Every excitation above 
this vacuum corresponds to a half-BPS state. Let the energies of the fermions in an 
excited state be Ei = eihw + where the Cj are unique non-negative integers. 

(The uniqueness is necessary because of the exclusion principle.) For the matrix 
model associated with the half-BPS states of A/" = 4 SYM we can simply choose 
u; = 1 in view of the conformal invariance of the theory. Each excited state can 
be parametrized by a set of integers representing the excitation energy of the z-th 
fermion over its energy in the ground state: 

r, = \{E,-Ef) = e,-i + l. (18) 

h 

The Tj form a non- decreasing set of integers r^r > tat-i > ■ ■ - Tj > 0, which can be 
encoded in a Young diagram T in which the z-th row has length r^. For example, 

1 

{4,3,1,1} 



(19) 

In our conventions the topmost row of the diagram is the longest. It is convenient to 
introduce another variable Cj with 1 < j < A^ so that [25] 

Cn = ri ; c^-i = r^+i - ; z = 1, 2, . . . , (A^ - 1) . (20) 

The variable Cj counts the number of columns of length j in the diagram associated 
to {r^, r^-i, . . . , Ti}. A useful relation is 



'i+i 



ei+i-i = CN-i + ■■■ + cn- (21) 



The fully anti-symmetrized wavefunction for the A^ fermion system can be expressed 
as a Slater determinant: 

■^Tiq) = -^det(vl/,^(g,)) ; vl/„(A) = A(n)i7„(A/v^) e-"'/2^ (22) 

Here i7„ is a Hermite polynomial, and A{n) normalizes the single particle wavefunc- 
tion. Happily, the state \E't associated to a Young diagram T turns out to be exactly 

^Strictly speaking we are describing an independent basis of operators in a U{N) theory. The 
vanishing trace of SU (N) matrices produces relations between some of the Or (see |2i]). At large- iV 
this will not change our results much. 
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the same state as created by the operator Or associated to the same Young diagram 
via ()17|) . This is because the asymptotic behavior of the Schur polynomial associ- 
ated to a Young diagram matches the asymptotic behavior of the Slater determinant 
13 II]. The conformal dimension of the operators, A, is equal to the number of boxes 
k appearing in the diagram T. 

Following the reasoning in Section 2, we are interested in the structure of typical 
half-BPS states of very large charge A = J = A^^. The energy in the excitation 
of the fermions is comparable to the energy in the Fermi sea. In analogy with the 
situation for Schwarzschild black holes, we will show that almost all such half-BPS 
states lie close to a certain limiting Young diagram. Because of the large amount of 
supersymmetry we do not expect these states to have a degeneracy large enough to 
be associated with a classical horizon in the dual spacetime description. Nevertheless, 
half-BPS states are related to extremal charged black holes that are on the verge of 
producing a horizon [S]. 



3.1 Structure of typical states 

Highly excited states of the large- free fermion system can reliably be studied in a 
canonical ensemble in which temperature rather than energy is held fixed. By stan- 
dard reasoning the error between the canonical and microcanonical approaches will 
be small in the thermodynamic (large- A^) limit in which we are principally interested. 
Thus we will use the free fermion formulation to study the typical shape of Young 
diagrams with A boxes. Related studies were carried out in jTH . 

Using the notation introduced at the beginning of Section 3, the canonical parti- 
tion function of fermions in a harmonic potential is given by 

Z = ^e-'^^= J2 e-'^^S»(^'+^/2) ^23) 

0<ei<e2---<ejv 

with u = 1 here. The restricted sums over the excitation energies Cj (or equivalently 
the row-lengths rj = Cj— 2 + 1) can be replaced by unrestricted sums over the variables 
Cj in (PU)) that count the number of columns of length j. Carrying out this trans- 
formation of variables and setting Z = Ze^^^'^"^ to remove the irrelevant vacuum 
energy, the partition function becomes 

oo _ N 

Z= Y: ^"'"^^^"^ = 1113^' (24) 

ci,C2,.--,cjv=l j=l 

with P = phw. Defining 

q = e-^ (25) 
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the temperature of the ensemble is fixed by requiring^ 

{E)=A = q^ log Z{q) = ^ . (26) 
^ j=i ^ 

In the large- limit we can approximate the sum on the right hand side by an integral, 
provided that the contribution from the limits is not too large. Assuming this to be 
true, as we will check self-consistently later, we have 

log(l-g^)-6Li2(g^) _ Li,(l-g^) 
6(logg)2 (logg)2 ^ > 

Here Li2(a;) = Xln^"/"^^ dilogarithm function. Since q = e~^, depending on 

how A scales with N, in the large- iV limit, q^ can approach 0, 1, or a finite number 
in between. 



A 



i V 



dj 



j=0 



jq^ 7i^ — 6N{logq) 
l-qi ~ 



On the scaling of temperature with A^: We can determine the scaling of the 
temperature in the ensemble with by estimating the number of microstates of a 
given conformal dimension A. The exact number of such microstates is given by 
the number of partitions of A into at most parts. (Equivalently, by the number 
of partitions of A into parts no bigger than A^.) We shall denote this number by 
p(A, A^), and so the entropy will be given by 

S = log p{A, N) . (28) 

Clearly, if A < A^, p(A, A^) = p(A) where p(A) stands for the unrestricted number 
of partitions of A. When A > A^, and in the large- A^ limit, we can still obtain an 
upper bound on the scaling of the number of states with A^ by applying the Hardy- 
Ramanujan formula [201 

^^logp(A) ^27rv/A76-logA-logy48 + C(l/VA). (29) 

Using this expression we can extract the leading A^ dependence in the temperature 
by computing the variation of the entropy with respect to the conformal dimension 

1 1 

^-^ + 0{A-'), (30) 

where we have neglected coefficients of order one, i.e. those not scaling with A^, 
since in any case the right partition problem of interest is not computed by p(A). 
Nevertheless this is a useful estimate because in large- A^, p{A) and p{A, N) will differ 

^Note that the temperature is simply a Lagrange multipher fixing the energy of our ensemble 
of half-BPS states. We are not studying a thermal ensemble in the Yang- Mills theory which would 
necessarily break supersymmetry. 
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by a factor of order unity. Thus we learn that for typical operators of conformal 
dimension A ~ A^^, the temperature should scale linearly with N. Equivalently, we 
are dealing with an ensemble in which approaches a finite number, even in the 
large- iV limit. 

Given the above microcanonical estimation, we shall assume that at large- iV we 
can write f3 = a/N for some constant a. Then 

A = N^ ( '^'^^',7""^ )^7iV^ (31) 

In other words, 

(3 = - =^ A = ^N\ (32) 

Having fixed f3 in this way we can compute the shape of the typical Young diagram 
with A boxes. 

The expected number of columns of length j is 

(c,) = ^ , (33) 
from which we derive the expected rowlength 

3=0 ^ 

In the large- limit overwhelmingly many half-BPS states will have associated dia- 
grams that lie arbitrarily close to the limit shape with rows of length (rj). Indeed, 
the diagrams that do not lie along the limit shape have vanishing likelihood. It is 
interesting that the columns of the Young diagram, which are interpreted as giant 
graviton D-branes [1^1 El 113 UHl IH CHI , are acting like a gas of bosons despite the fact 
that we are studying the excitations of a fermionic system. It is well-known that small 
fluctuations of a Fermi sea in the c = 1 and AdS/CFT contexts can be bosonized 
to described small bosonic fluctuations of spacetime [131 HIl I2H] • We are explaining 
an interesting kind of bosonization that produces bosonic D-branes from collective 
fermionic excitations. 



Entropy: The entropy of the above canonical ensemble can be obtained in the usual 
way by means of a Legendre transform of the free energy. This is a straightforward 
calculation, and the result is that for the scaling regime with A = 'jN^ and 7 kept 
fixed as — > 00 the entropy becomes 

S{A, N) = -l^-- U,{e--)j + a-, (35) 

with a the solution of ()31|). To get an idea of the extent to which the constraint 
on the number of rows reduces the number of partitions, we take the example with 
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7 = 1, i-e. A = N"^, and compare S{A,N) with the entropy S'o(A) = logj9(A) of 
unconstrained partitions with A boxes. For 7 = 1 one can numerically solve (|HT|) and 
one obtains a ~ 0.81465. Inserting this into (j35p yields S ~ 2.214iV. The number of 
unconstrained partitions is given by the Hardy-Ramanujan formula 



e^^^ / 1 \ vr 

p{A) ^ 1 ^ ; c = ^. (36) 

4v^A V Ic^fA) V6 

Thus So (A) ~ 27rA/A/6 ^ 2.565A^, which shows that fixing the number of rows to 
be iV cuts down the entropy by a factor of about 0.86. In any case, this numerical 
analysis shows that indeed the dominant dependence in the entropy is correctly 
captured by the Hardy-Ramanujan formula as above. 

As an aside, we briefly compare this to the asymptotic formula 



p(A, N) = p{A) exp ( -^e-"^ ) (37) 



where this is valid in the limit where A^, A 00 keeping 



A^ logA 



fixed. ^ In this regime A^ grows faster than a/A, namely A^ grows as vAlog A instead. 
If we nevertheless insert A = A^^ in ()37|) we obtain an entropy S ^ 2.349A^ which 
is 0.92 times the unconstrained entropy. Either way, we find that the numbers of 
constrained and unconstrained partitions are quite close to each other. 

Temperature: We can also formally derive a temperature by taking the derivative 
of entropy with respect to the energy A. Since q = exp{—a/N), it comes as no 
surprise that this yields 

T~ dA~ N ^^^^ 
There is no physical temperature here as we are dealing with half-BPS states. Rather 
T should be understood as the effective temperature of a canonical ensemble that 
sums over half-BPS states with energy sharply peaked at A}^ 

Limit shape: Let us introduce two coordinates x and y along the rows and 
columns of the Young diagram. In our conventions, the origin (0, 0) is the bottom 
left corner of the diagram, and x increases going up while y increases to the right. 



^By keeping x fixed here we mean that this parameter does not scale with TV as the limit is taken, 
but is otherwise free to rmi from —00 to 00. 
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One can show more generally that if A scales as N^, the temperature goes as N^/'^ . 
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In the fermion language, x labels the particle number and y its excitation above the 
vacuum. The expression (j?H) can then be written as 



N 

y{x)= J2 • (40) 



i=N-x 

Taking the expectation value, 

N 



(y(x)) = y{x) = (c^) • (41) 

i=N-x 

In the large- limit, we can treat x and y as continuum variables, and the summation 
becomes an integral^^ 

pN 

y{x) = / di (q) 

J N-x 

log(l - e-^^) log(l - e-^^^-"^) 

= C(AiV)-i^SU_^ 1 ^42) 

After some elementary algebra, this curve can be rewritten as 

g^-^ + t{p, N)qy = 1 ; t{(3, N) = . (43) 

The precise meaning of the limit shape is that in the large- limit, and after 
a suitable rescaling of the Young diagram, the boundary of the Young diagram ap- 
proaches the limit curve with probability one. In other words, in the large- iV limit 
almost all operators with conformal dimension A will have associated Young diagrams 
that are vanishingly small fluctuations about this curve. Limit shapes for Young di- 
agrams have also been discussed by Vershik as part of an extensive mathematical 
literature on random partitions of integers. 

When we approximate the sum in ()42|) by an integral, one may worry that we 
miss some important aspect of the physics. For example, it could happen that {cn) 



-'^-'^More carefully, we take N oo while /i — > such that the Fermi level HN is kept fixed. 
Rescaling x and y appropriately leads to the integral. 

^^Our limit shape H43|l differs from the limit shapes of Young diagrams that appear in explicit 
instanton calculation in A/" = 2 super Yang-Mills theory and in computations of certain Gromov- 
Witten invariants PJj- This is because we are averaging over diagrams using the uniform measure 
rather than the Plancherel measure. We are treating all the partitions as a priori equiprobable, 
rather than assigning weights to partitions based on the dimension of the representation of the 
associated Young diagram. This is the correct procedure in the present case because the volume 
of the residual orbits of the symmetric group Sn that remain after gauge fixing the original U{N) 
symmetry have been accounted for in the way we parametrize free fermion excitations. 
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becomes macroscopically large of order VA (say, if for instance there were some 
version of Bose- Einstein condensation in the theory). For our ensembles, this only 
happens in certain limiting cases, which we deal with as necessary. 
It is convenient to rewrite the limit curve as 

q-q^ + qy (l - g^) = 1 . (44) 

As we have seen, if A = 7A^^, the inverse temperature goes as (3 = a/N. The 
above curve then becomes: 

e-"g-^ + (l-e-")g^ = 1 . (45) 

The limit shape simplifies in two natural regimes. First, we focus on the base of the 
limit shape, that is, at x ^ A^. In this regime, the limit curve behaves like a straight 
line: 

— Ct 

y{x) = X, |r < 1 . (46) 

Secondly, we focus on the top of the limit diagram. Introducing the variable x = 
N — X 'ti N, we find logarithmic behavior: 

Fluctuations: The accuracy of the limit shape as a description of typical states 
is determined by the size of fluctuations in individual states. The variance in the 
number of columns of length j turns out to be 

Var(c,) = ((c, - (c,))^) = {c'^) - {c.f = = ^ (48) 

just like the fluctuations of population numbers in a bosonic gas. Thus the standard 
deviation divided by the mean is 

(49) 





f 1 ^ 




f 1 ^ 


1 ' 


(9) ' 











where the last equality sets /3 = a/N. Thus, we are seeing the standard -y/n random 
fluctuations in counting populations, enhanced by a factor representing the tendency 
of bosons to condense to low-energy states. This factor matters most for the number 
of short columns (small j). 

It is also interesting to compute the fluctuations in the limit shape ()44|1 . In order 
to calculate this using ()4H). we use the fact that the numbers of columns of different 
lengths are uncorrelated variables: 

{c^Cj) = {Ci){Cj) i^j. (50) 
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Using this in the definition of the variance gives 

Var?/(x) = ((y(x) - {yix))^) = ^ Var(Q) 



N 



(51) 



i=N-x 



Taking the continuum hmit as in (j42j) we find 

1 1 



Var 



+ 



-l3{N-x) ■ 



(52) 



We can assemble this into a probabihty distribution describing the hkehhood of small 
fiuctuations around the limit shape 



Pr[y(x)] = C exp 



N 



dx 



y{x)f 



2Var y(x) 



(53) 



The normalization constant C is a functional determinant of Var?/(x) in the usual way 
and ensures that J T>y{x) Pr[y(x)] = 1. After a suitable conversion of the variable 
y(x) into a probability distribution, this expression could be compared to Sanov's 
theorem Q which described the probability of having a large deviation for a particular 
Schwarzschild microstate from the typical state. 



3.2 Error estimates 

It is useful to estimate the error incurred in our considerations by replacing sums by 
integrals and by taking a canonical rather than a microcanonical approach. Different 
ensembles typically agree at leading order in an expansion in 1/E, but the subleading 
terms will typically be different. To get an idea of the magnitude of these errors, 
we turn to the simple case of diagrams with A boxes and no constraints on the 
number of columns. This is reasonable to study because, as we showed in Section 3.1, 
the restriction on the number of rows makes only a minor difference to the overall 
degeneracy of states. 

Taking the logarithm of the Hardy-Ramanujan formula ()36|) . we see that 

logp(A) = 27rv/A76-logA-logy48 + C(l/VA). (54) 

When we replace sums by integrals, we recover only the first term in this expansion, 
and from this we see that the typical error introduced by replacing sums by integrals 
goes as log E/\fE. Below we will see another example where the difference between 
the sum and integral yields an error of order 1 / \fE. Though we have not extensively 
studied the errors in all ensembles discussed in this paper, it seems very likely that 
similar estimates apply there as well. 

In a canonical analysis, one computes the partition function Z{q) = ^^]9(A)g^: 

Z = ^ exp (^A logg + 2n^/A/6 - log A - log ^48 + 0(1/ VA)) . (55) 

A 
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This partition function peaks at A = Aq, where 



logg = --J= + + OiA-'^'). (56) 



^/6A^ A, 







The variance is obtained by differentiating logZ twice at A = Aq and then taking 
the inverse: 

logZ 



9A2 



(57) 



A=Ao 



This works out to be 



^,^y24A|^24A,^ (58) 

IT TT^ 

As a check, one can compute o"^ simply via the standard expression for the variance 
of the energy in a canonical ensemble 



2 



a 



n^q"" 2C(2) 



n>0 



which reproduces the leading term in (j^H|) when we use (jSSI)- From the above we infer 
that in this case the mistakes that are made when replacing sums by integrals are all 
of order 1 / \fK. 

A slightly more complicated example is the expectation value of the number of 
columns with j boxes. The canonical ensemble answer is very simple, 

(c,> = (60) 
The micro canonical answer is on the other hand given by^^ 

To evaluate this, we use the asymptotic form of p(A) in (jHH), and expand the exponent 
in kj. We discover that 

5^p(A - kj) ~ J^exp (-vr^ + X + ' ' ') ' ^^^^ 



A:>0 fc>0 

where the dots indicate terms of higher order in kj or terms suppressed by higher 
powers of A. If we keep only the leading term, we can do the sum exactly and recover 
fjHnjl with q = exp(— tt/V 6A), the leading answer to (jKHjl . The two obvious sources 
of mistakes are dropping other terms linear in kj, which will shift the value of g to a 

""^^A given partition with precisely m columns of length j will appear a total of m times in the sum 
evaluated in the numerator, as required. 
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slightly different value, and dropping the quadratic and higher terms in kj. It seems 
that the main source of error arises from the interference of quadratic terms with 
linear ones. But these barely have any impact unless j starts to become of order 
A^/^, so that for A; = 1 there is some suppression from the quadratic terms. Thus we 
reach the preliminary conclusion that the microcanonical answer and the canonical 
answer differ by a factor of the form 

^+°(a^-;^0- '''' 

We have checked numerically that this is the right behavior. 

The next step is to consider the variance in the distribution of the number of 
columns of a given length. In Section 3.2 we showed that this was 



(1 - q^r 

The exact microcanonical answer is^^ 



- TT^.- (64) 



,2. _ E.>o(2fc - 1)P(A - kj) 

(S) - ■ (65j 



With the same manipulations as above we find that this equals 



q-' {! + (/) , j - ^^ 



so finally we find that the microcanonical answer for o"^ is given by 



The corrections are found in the same way as above. The spread in the expectation 
value of Cj is therefore 

1-'^'(^ + 0{-^,,4^)). (68) 



(9) ^ V ^A3/4'yA' 

This illustrates how the number of columns of each length is distributed as we probe 
all the states in a microcanonical ensemble. 

In the subsequent considerations in this paper we will ignore corrections such as 
the ones studied in this section, but it is important to keep in mind that they are 
there. 



^^A given partition with m columns of length j will appear in the sum in the numerator once for 
each 1 < k < m. Thus it will contribute a weight X^^i(2fc — 1) = m^, exactly as required. 
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3.3 The superstar ensemble 

In this section, we analyze an ensemble of Young diagrams in which not only the 
number of boxes A and the number of rows are held fixed, but also the number 
of columns Nq- Physically, the last constraint amounts to holding fixed the number 
of D-branes in spacetime ^21 I3|-^^ One might expect that this ensemble includes 
a microscopic description of the extremal "superstar" configuration since this 
spacetime involves a constraint relating the total energy, the amount of flux and the 
number of D-branes. We will indeed find a description of the superstar once we 
restrict the three parameters to obey A = NNc/2. For arbitrary values of the three 
parameters we find a finite "temperature" deformation of the superstar. 

The addition of the column number constraint changes the statistical properties 
of the system considerably. In the previous ensemble, the average energy could grow 
without bound by increasing the temperature of the ensemble. The larger the energy, 
the larger the entropy of the system. In the present set-up, this is no longer true. 
Given an ensemble characterized by the pair {N, Nc), the conformal dimension A 
lies in a finite interval : 

{E) = Ae [Nc, N-Nc]. (69) 

The lower bound corresponds to a one row Young diagram with Nc boxes; whereas 
the upper bound corresponds to a rectangular diagram in which all rows have Nc 
boxes. Clearly, both bounds have a unique microstate, and so we can conclude our 
system has vanishing entropy in both situations. 

Intuitively, one expects that as one increases the energy slightly above Nc, the 
entropy increases, heating up the system to a finite positive temperature. This be- 
havior should continue until the entropy reaches a maximum, which is achieved at 
As = Nc{N + l)/2, when half of the allowed boxes are filled in. If the energy goes 
beyond this value, we expect the entropy to start decreasing as the degeneracy of mi- 
crostates will decrease. Indeed, there is a one-to-one map between Young diagrams 
with A boxes and Nc{N + 1) - A boxes, and therefore 5'(A) = S{Nc{N + 1) - A). 
This shows that in this regime our system will have negative temperature.^^ As the 
energy continues to grow, the absolute value of the temperature will decrease, reach- 
ing a zero value for the upper bound. This discussion strongly suggests that the 
system achieves an infinite temperature when the entropy is maximized. Below we 
will see that this is precisely what happens, and in Section 4 we will show that the 
maximum entropy configuration is described in gravity as the superstar of [S]. 

The description of typical states will be based on the canonical ensemble analysis 

^^At a more technical level, introducing the parameter Nc allows us to have some control over 
the way in which this number scales with N. 

i6\Ye recall that the temperature in these ensembles is not physical. It is merely a parameter to 
characterize typical states, which in this case allows us to correctly describe the situation in which 
an increase in energy causes the entropy to decrease. 
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in Section 3.2, with an additional Lagrange multiplier to implement the constraint 
on the number of columns in the Young diagram. Using the same notation as before, 
the partition function for the system, once we have removed the energy contribution 
from the vacuum, is: 



oo 



J2 e-''^^^'^~^(^^'^-''^) , (70) 



ci,C2,---,cjv=l 



where A stands for the Lagrange multiplier. (Recall that Cj counts the number of 
columns of length j.) Thus we will require the condition 

— = ^ J]c, = iV^. (71) 
j 

Proceeding as before, this partition function can be written as 

j 

The temperature of the ensemble is fixed by requiring 



N 



{E) = A = qd, log Z{C, q) = J2 {zfl ' (^3) 
whereas the Lagrange multiplier (A), or equivalently C, will be fixed by 

N N ^ j 

where we already computed the expected number of columns of length j, i.e. {cj). 

We shall examine the different natural temperature regimes in this ensemble, to 
check the claims made in our preliminary discussion. For the time being, we will 
keep fixed and vary j3. We can also take scaling limits involving both (3 and A^. 
These two ways of taking limits are not identical. They lead to similar qualitative 
conclusions but the quantitative details are different. 

So for now we keep N fixed and consider the large-fi expansion. If the temperature 
is positive., this corresponds to the approximation q = e~^ -C 1. In such regime, both 
sums in (fTSj) and (fTIj) are dominated by their first term. This allows us to conclude 

A ~ iVc ~ (ci) = , (75) 

1 - Cq 

which indeed corresponds to a diagram with a single row of Nc boxes since the 
dominant contribution comes from columns of length one. 
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If the temperature is negative, corresponding to the approximation q = e ^ ^ 1, 
the number of columns 

Nc - {cn) = , (76) 

is dominated by the number of columns of length A^. This is so because the function 
x/{l ~ x) is monotonically increasing in x. Thus the corresponding diagram is a 
rectangular one, as expected. To further check this interpretation, we can compute 
the dominant contribution to the energy ()73p 

oo N 
s=l j=l 

oo 

~ iV^(Cg^r = iViVc, (77) 



which indeed agrees with the energy of a rectangular diagram with rows of Nc 
boxes each. 

Let us analyze the more interesting high temperature regime, or equivalently, 
the small-/? expansion. We shall treat both the positive and negative temperatures 
together including a sign in the definition of the expansion parameter (jS). We are 
interested in computing the entropy to second order in j3, to provide evidence for the 
existence of a maximum at /3 = 0. Using the identity, 

TV 

S = (3A + \ogZ = f3A- iVclogC - $^log(l -(q^) , (78) 

i=i 

we realize that we need to work out the expansions for ( and A at second and first 
orders, respectively. Expanding (f?^ and inverting the corresponding equation, we 
can find the expression for (: 



CiN, Nc) = -^[l+f3j;^+ P'^^iiN + 2) - u{N - 1))J + 0{(3') . (79) 
In the above expression we have introduced the notation: 

Nc , N(N + 1) , , 

Expanding (|7H)l . the conformal dimension at first order is 

A = ^Ai-/3 4^a;(l + ^)(iV-l) + 0(/?2) . (81) 
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Notice that at large- A^, the dominant energy is the one computed for the superstar 
in [5], suggesting that in the T — > cx) hmit the large- typical state is described in 
spacetime as the superstar. We will show this explicitly in Section 4.3. The typical 
states in this infinite temperature ensemble correspond to Young diagrams that are 
nearly triangular. That is, on average, there is a constant gap {uj) of energy between 
the excitations of the (j + l)-st and j-th fermions, the first fermion having an average 
energy uj itself. Moreover, having computed the linear j3 dependence allows us to 
confirm that for positive temperatures, the conformal dimension is smaller than the 
superstar energy, whereas it is for negative temperatures that the conformal dimension 
can exceed the latter, in agreement with our general entropic arguments. 

Finally, if the superstar ensemble maximizes the entropy, the entropy expansion 
in f3 should have no linear dependence in it and its second order coefficient must be 
negative for all values of A^, Nc- Carrying out the computation, we obtain 

^ - - ixi^S^) - - 1) + ■ 

Since > 1, the coefficient of is negative, for any value of (A^, Nc), as expected. 
According to our discussion, we should identify the first term above as a microscopic 
derivation for the entropy of the superstar. 

Limit shape: Proceeding as in the previous section, in the large-A^ limit, we can 
deduce the continuum limit curve: 

N „N 
j=N-x 

= ■^iog(i-Cg"')--^iog(i-Cg"^-1 

= c)--^iog(i-Cg^-1 , (83) 

where our conventions for both {x, y} along the Young diagram and its origin are as 
before. This limit shape can be written more symmetrically as 

^ qN-. + ^ ^y-Nc ^ 1 ^ ^ ^ qNa-C(/3, () ^^^^ 

Since we are working in the limit of large A^, Nc the slight asymmetry between A^ 
and Nc that we had before is irrelevant. In particular, the points (x, y) = (0, 0) and 
(x, y) = {N, Nc) lie on the limit curve. Therefore 

' ' C = -^^- (85) 



qN _ q-Nc ' ^ qN _ q~Nc ' 

This is equivalent to the constraint defining the ensemble (f?^ as one may easily 
verify. Exactly the same limit curve also appears in theorem 4.7 in [22], where it 
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is shown that this is indeed the hmit curve if we take A ^ oo keeping N/^/A and 
Nc/^/A fixed. 

To make the symmetry between N and Nc more manifest we rewrite the hmit 
curve (jS^ as 

ag^-^ + (3qy = 1, (86) 

with 

l-q^c 

" = ^ = 1 _ qN+Na ' ^ = 1 _ qN+Na ' (§7) 

The energy, or in other words A, can be computed from ()73|) which is of course the 
same as computing the area under the hmit curve. We obtain 

Li2(a) + Li2(/3) + logalog/3-7rV6 
(logg)^ 

which is indeed symmetric under N ^ Nc- The parameter q in the grand canonical 
partition function is nothing but e~^^'^ and we saw above that as T — oo the entropy 
reaches a maximum and is dominated by triangular Young diagrams. It is easy to 
confirm this from the limit shape: expanding eq. ()86|) around g = 1, we find that the 
limit shape becomes 

^|^(,-l) + O((,-lO = 0. (89) 

To leading order this is indeed a straight line 

Nc 

y = -j^x = ujx (90) 

whose slope is ratio of the number of columns to the number of rows, similarly to 
what we saw above. 

The Young diagram that is encoded by the limit shape (jUUj) admits two natural 
interpretations. We can read the diagram from left to right or from top to bottom. 
If we read from left to right, the columns represent giant gravitons, D3-branes that 
wrap C 5*^. The maximum angular momentum of any of these giant gravitons is 
Nr < N but their number is unconstrained. The geometry given by this limit shape 
arises as a bound state of Nc giant gravitons. The angular momentum of the y-th 
giant graviton is A^^^; — yNn/Nc- If we read from top to bottom, the rows represent 
dual giant gravitons, D3-branes that wrap C AdSs. We can have at most Nr < N 
of these because dual giant gravitons are stabilized by fiux, and the fiux through 
the last dual giant represented by the bottom row of the Young diagram is iV — Nr. 
The same geometry will bound state of the Nr dual giants. The angular 
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momentum of the x-th dual giant graviton isy = xNc /Nr}'^ As we shall see explicitly 
in the next section, this geometry is the superstar. We defer further discussion of this 
for the moment in order to further explore the thermodynamic properties of this limit 
shape. 

We briefly return to the entropy of the system. The scaling limit relevant for the 
limit curve is one where A ^ oo while keeping A = \/~K/Nc and yU = \/~K/N fixed. 
In this regime q will behave as exp(^/A/A) for some fixed ^, and a and (3 as defined 
in (jS7|) will be finite as well. Given A,yU, the three functions are the solutions 

of the three equations 



a 

e 



1 _ eC^+A')? 

1 -e^« 
1 _ e(^+/^)5 

Li2(a) + Li2(/3) + logalog/5 - 7rV6. 



(91) 



and the entropy that we obtain by Legendre transforming the free energy is 

^(A,/i,A) ~ (-2^- Aloga-/ilog/5)yA (92) 

These equations are difficult to solve in general. However, the linear limit shape is 
easy to get by taking ^ ^ 0, and we see that a = A/(A + yu) and /5 = yu/(A + /x). For 
the linear limit shape with Nq columns and rows we therefore find that 



A = ^ATATc, 



A 



c 



N 



2N 



(93) 



and the entropy is given by 

S ~ log 

which agrees with the leading term in 



(A + /i) 



A^/i^- 



(94) 



4 1/2-BPS States: Effective description in gravity 

In the most general type IIB geometry invariant under 5*0(4) x 50(4) x M 
preserving one-half of the supersymmetries, consistent with A = J, is determined. 

^^That there are two such interpretations of the same background is a manifestation of the parti- 
cle/hole duahty of the quantum HaU description of the half-BPS states [3311^. In fact this exchange 
symmetry even manifests itself away from the BPS hmit. It is shown in |85j that the entropy of the 
non-extremal superstars with charge G in a background with N units of flux is equal to the entropy 
of charge N superstars in a space with G units of flux. This is a non-extremal incarnation of the 
exchange symmetry between rows and columns in the ensembles we are considering. 
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We develop the correspondence between these geometries and effective semiclassical 
description of half-BPS states in the dual Yang-Mills theory, particularly those space- 
time(s) that describe states with A ~ iV^. We argue that almost all such states have 
as an underlying structure a "quantum foam" whose universal effective description 
in supergravity is a certain singular spacetime that we shall dub the "hyperstar". 
The singularity arises because the classical description integrates out the microscopic 
details of the quantum mechanical wavefunction. 

In all the relevant half-BPS supergravity backgrounds are constructed in terms 
of metric and self-dual five-form field strength F(5) = F A dfl^ + F A dQ^, where dQ^ 
and dQs are the volume forms of the two three-spheres where the two 5*0 (4) s are 
realized. The metric is: 

ds'^ = -h~^{dt + Vidxy + h'^{dri'^ + dx'dx')+T]e^dnl + 7]e-^dnl (95) 
h-^ = 2ricoshG, (96) 
T]drjVi = eijdjZ, vidiVj - djVi) = ei^d^z (97) 

z = ^tanhG (98) 

where z = 1,2 and *3 is the fiat space epsilon symbol in the directions ?7,xi,X2. In 
addition there is a self-dual flux that depends on the function z. The entire solution 
is evidently determined by the single function z, which obeys the linear differential 
equation 

dAz + vdrj(^^^ = 0. (99) 

The function $(//; x^, x^) = zr]~'^ satisfies the Laplace equation for an electrostatic 
potential in six dimensions that is spherically symmetric in four of these directions. 
The coordinates Xi,X2 parametrize an M^, while rj is the radial coordinate in the 
transverse to this auxiliary six- dimensional manifold. Note that xi and X2 have 
dimensions of (length)^. The analogy with electrostatics allowed to solve for the 
bulk geometry in terms of the boundary condition z{0, x[, at the origin (77 = 0): 



z{v; xi, X2) = — dx[ dx'2 Yi '-w^. — ^ ■ (100) 

z(0; a;i, X2) = ±1/2 are the only boundary conditions compatible with non-singular 
geometries. In terms of the asymptotic AdSs x geometry, at points on the (a;i, X2) 
plane where z = 1/2 (—1/2), an in the (AdSs) shrinks to zero size. This 
topologically complex two-plane is spliced smoothly into the full bulk spacetime. Of 
course, other boundary conditions are also compatible with a half-BPS condition; 
these will give rise to singular spacetime geometries. Notably, |2(0;xi,X2)| > 1/2 
leads to closed timelike curves while \z{0; Xi,X2)\ < 1/2 gives null singularities [36, 37J. 

Following jni; in order to match these solutions with states in the field theory, 
consider geometries for which the regions in the {xi,X2) plane where z = —1/2 are 
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compact. We will call these regions collectively the droplet V. The quantization of 
flux in the geometry leads, even in the semiclassical limit, to an identification between 
h in the dual field theory and the ten-dimensional Planck length: 

h ^ 27r£j. (101) 

(The units in this identification are unusual because of the dimensions of {xi,X2)-) 
Furthermore, the area of the droplet T> is quantized and must equal A^, which is 
equivalently the total amount of five-form flux, the rank of the dual gauge group, and 
the number of fermions in the free Fermi picture of half-BPS states: 

^^-IM' ''''' 

where we used ()101|) . The conformal dimension A of a given configuration (state) is 
computed by 

d X 1 3/ -j^ ~|~ ^2 ^ f f d /'-ir\o\ 

(103) 



27ih2 h 2 Vie 27rn 
These equations have a remarkable interpretation in terms of the hydrodynamic limit 
of the phase space of the fermionic system used to construct half-BPS states in Section 
3 jHllll. Specifically, if we identify the {xi,X2) plane with the phase plane of a single 
fermion in a harmonic potential, 

{xi,X2) ^ (g,p), (104) 



then in the hydrodynamic limit for an A^-particle fermionic system ()102|1 and ((n 
precisely compute the number of fermions and the total energy over the vacuum. 
In fact, the first term in ()1()3|1 calculates the total energy of a fermion droplet in 
phase space while the second term in ()103|) subtracts off the vacuum energy. This 
identification gains support because the droplet corresponding to the filled Fermi sea 
{i.e. a circular disk) reproduces the geometry of the AdSs x vacuum. Based on 
such observations [13] proposed that the {xi,X2) plane of their half-BPS geometries 
should be identified with the single particle phase space plane of the dual fermionic 
system. 

For our purposes it is convenient to introduce the function 



m(0; xi, X2) = ^ - ^(0; xi, X2) , (105) 

which takes values one or zero for non-singular geometries. For these, the energy 
integral ()103|1 may now be written as 



N = I ^u{0;xuX2). (107) 
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These expressions resemble expectation values, suggesting that the function u should 
be identified with the semiclassical limit of the quantum single-particle phase space 
distributions of the dual fermions. (Also see jSHl ISHl HOI which appeared while this 
paper was being prepared.) We will see that u does not always assume the values 
0, 1; the generic spacetime description of half-BPS states will be singular. 

4.1 Proposal 

Below we shall carry out a detailed study of the semiclassical limit of the phase 
distributions of the half-BPS fermionic system and explain how the function u arises as 
a coarse-grained effective description of underlying quantum dynamics. However, even 
before doing this analysis, it is possible to propose a detailed map between specific 
Young diagrams representing half-BPS states and the corresponding semiclassical 
geometries, which will all be circularly symmetric in the Xi, X2 plane. Consider states 
whose descriptions in terms of Young diagrams have a well-defined semiclassical limit. 
For the moment we will take this to mean that as ^ — with A^^ fixed, the perimeter 
of the Young diagram approaches a smooth curve y{x) with x parametrizing the rows 
and y the columns. Examples of such curves were derived in the previous section 
as limiting shapes of typical partitions in various thermodynamic ensembles. We 
propose that in such situations the integral formulae (|1(J6|) . (I1U7|) extend to differential 
relations: 

dx , (108) 

{y{x) + x)dx . (109) 

Here we have written m(0; Xi, X2) as ^(0; r^) in terms of the radial coordinate r in the 
{xi,X2) plane in view of the U{1) symmetry present for single Young diagram states. 
In terms of the phase space interpretation of ()106|) . ()107|) . the first equation simply 
relates the number of particles in phase space within a band between r and r + rfr to 
the number of particles as determined by the rows of the associated Young diagram. 
The second equation matches the energy of the particles in phase space within a ring 
of width dr to the energy in terms of the Young diagram coordinates. 
Combining equations ()108|) and ()109p . we find that 

y{x) + X = r"^ / {2n) . (110) 

This is just a check on the consistency of our proposal, since indeed y{x) + x is the 
total energy at a given point x in the diagram, and so must match the contribution to 

^^Recall that xth fermion has an excitation energy of hy(x) over its ground state energy hx. 
Note that dx hx — , which is indeed the total vacuum energy of N fermions in a harmonic 
potential. Equivalently, x dx is the continuum version of the vacuum energy stored in fermion x. 



u(0; 



dr' 



-dr' 
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the energy r^/ {2h) at a given radial distance r in the (xi, X2) plane. Taking derivatives 
with respect to x and using (jlOSj) . we derive the relation 



1 + y' 



1 



^ z(0;r^) 



l y'-i 

2y' + l 



(111) 



The right hand sides of the equations ()111|) should be understood to be functions 
of obtained by substituting the known boundary curve y{x) into (jllOj) and then 
inverting to find x as a function of r^. We thus establish a dictionary between the 
functions u{0; Xi, X2) that completely determine classical half-BPS solutions and the 
slope {y'{x) = dy/dx) of the Young diagram of the corresponding field theory state. 
The only ingredients we have used in making this identification are: (a) a differential 
relation between the number of flux quanta in gravity and the number of fermions 
in the dual fleld theory, and (b) a differential relation between energies on the two 
sides. Although it is suggestive that (xi,X2) look like phase space variables and in 
particular looks like the Hamiltonian, we did not use this in any significant way. 
Remarkably, in Section 4.2 we will show that in the semiclassical limit, at scales large 
compared to the exact quantum mechanical phase space distribution for the half- 
BPS fermion system precisely reproduces (llllj) . The distance from the origin is an 
energy. 

Despite the large amount of evidence we will present for our proposal, one could 
imagine going further and trying to prove it directly in the context of AdS/CFT 
correspondence. That the {xi,X2) plane can be identified with the phase space of a 
fermionic system was already demonstrated in [351 HUl 1^ • From this one could go on 
and try to show that the energy-momentum tensor of the dual matrix model, which 
couples to the bulk metric, indeed effectively integrates out all but one of the fermions, 
and one could perhaps show that higher, massive string modes are necessary in order 
to reconstruct the detailed structure of an arbitrary iV-fermion density matrix. This 
would emphasize once more that gravity is just an effective IR description of the 
microphysics, and that it is pointless to take the metric seriously once the Planck 
scale is reached. We leave such explorations to future work. 

4.2 Phase space distributions 

Above, we have discussed the analogy between (jl(J6|) . (jl(J7p and the phase space for- 
mulae for energies and populations of fermions in a harmonic potential. In this section 
we shall explore in detail the exact quantum mechanical phase space description of 
half-BPS states in Yang-Mills theory. In Section 4.3 we will extract the semiclassi- 
cal limit that is appropriate for comparing with the geometries of jT^j, but the exact 
quantum structures described here will contain lessons about the origin of singularities 
in classical gravity. 
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4.2.1 Wigner Distribution 

Given an n- dimensional quantum mechanical system described by the density matrix 
the Wigner distribution function on phase space is defined by 

1 /"°° 

W{q;p) = -—- / rfy (g-y|p„|g + y}e2^^^-^^/^ (112) 

This function is real and its projections give probability densities on the configuration 
and momentum space: 

j dqW{lp) = mn\p) , (113) 

df>W{q,p) = {^pM , (114) 
J dq J dpW{q,p) = tr(p„) = n . (US) 

Despite this, and except in special cases like Gaussian wavepackets, the Wigner func- 
tion itself is not positive definite and is thus not in general a joint probability distri- 
bution on phase space. Nevertheless, the Wigner function computes operator expec- 
tation values through the relation 

tT{pnAwCQ,P)) = j dqdp A{q,p) W{q,p) , (116) 

where the operator Awiq,f)) is Weyl (symmetric) ordered in q and p. (Equivalently, 
for particles in a harmonic potential it is symmetric ordered in the ladder operators.) 
The absolute value |H^(^,p)| is bounded. In general the Wigner function oscillates 
rapidly inside the classical torus and decays exponentially outside of it. Useful reviews 
include jl21 iSl El • 

As we discussed in Section 2, the half-BPS states in SU (N) Yang-Mills theory are 
related to the excited states of fermions in a harmonic potential. The corresponding 
phase space is 2A^-dimensional. However, the data relevant for the half-BPS solutions 
in gravity appear to involve only a two-dimensional phase space containing an effective 
single-particle distribution. One computes such an effective one-particle distribution 
for an A^-particle system by integrating out all but one of the particles from the 
A^-particle density matrix: 

/oo 
dq2dq3 ■ ■ ■ dqn {q2, ■ ■ ■,qn\pn\q2, ...,qn)= tr-?^2®. .. ®W„ (Pn) • (117) 
■oo 

Though it may appear that we unnaturally singled out one of the fermions in the 
above expression, due to the antisymmetry of the iV-fermion wavefunction pi is in 
fact completely symmetric in all the fermions. We are interested here in the Wigner 
distribution of a half-BPS state encoded by a single Young diagram with row lengths 
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TZ = {ri, ■ ■ - tn}- In the fermionic representation of the state, the excitation levels 
are JF = {/^ = n, ■ ■ ■ /^v = r^r + iV — 1}.^^ Each fermion is in an orthonormal pure 
state (\E'/.(g)) and the wavefunction of the iV-fermion system is given by the Slater 
determinant of the single particle wavefunctions 

^(^) = ^det(%(g,)) . (118) 

Due to the orthogonality of the the single particle density matrix becomes 

p^ix,y) = j dq^*{x,q)^{y,q) = \C\' V^}(x) z^/(i/) j , (119) 

where |Cp is a normalization constant. The effective single particle Wigner distribu- 
tion function is defined as 

1 r°° 

W{q,p) = -- / dy {q - y\p,\q + y)e''Py/'^ (120) 

In our case, we are dealing with fermions in a harmonic potential. The single 
particle wavefunctions are 

^/(g) = ^(/)^^/(g/v^)e-^'/^^ (121) 

A{f) is a normalization factor and Hf is a Hermite polynomial. Following our general 
discussion, and using the identity 

/oo 
dxe-^' Hmix + y)H4x + z)=2''y/^m\z''-'^L';;;"^{-2yz), m < n , (122) 
■oo 

the single particle Wigner distribution turns out to be 

W{q,p) = ^e-(^^+^^)/^5^(-l)^L, {liq'+p')) , (123) 

where Ln{x) is a Laguerre polynomial (the normalization is L„(0) = 1). The distri- 
bution is invariant under rotations in the two-dimensional phase space as it should 
be for the states under consideration. 



A single excitation: Since the final distribution function is the sum of single 
particle Wigner distributions, we begin by considering the Wigner distribution of 
single fermion at excitation level n |45j: 

WrM.P) = \{-ire-'''I^K{AH/h), H{q,p) = hp' + q') . (124) 
^^We have dropped the irrelevant h/2 zero-point energies. 
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Figure 1: The Wigner distribution function of a single fermion with excitation level 
n = 50 as a function of H/ h. 

At the origin of phase space, W„(0) = {—!)"' /nh. As we move in the radial direc- 
tion, the Wigner function oscillates with increasingly broad peaks (see Fig. . The 
positions of the maxima and minima are given by the solutions to the equation 

(^-| + L„(C) - nL„_i(C) = 0, (125) 

and thus the oscillations occur at a scale \/h. The final extremum is always a max- 
imum and is positioned at C ~ An, which coincides exactly with the classical orbit. 
Beyond this maximum, the function decays exponentially. Thus, the Wigner distri- 
bution is sensitive to fine quantum mechanical structures at the h scale. In particular, 
the Wigner distribution for a single fermion is not well localized at scales below the 
classical energy.^° Indeed, it is not even positive. 

Although the Wigner distribution is the most natural distribution in phase space 
the fineness of the quantum mechanical structure and the lack of positivity 
mean that it is not a good candidate for mapping into the boundary data required 
to determine the gravitational dual. In particular to obtain the classical geometries 
described in |1B], we must ask how to take the appropriate semiclassical limit of the 
Wigner distribution. The semiclassical limit requires — > 0. To keep the Fermi level 
fixed we should take ^ oo in such a way that hN remains fixed. Since energy levels 
are spaced by h, phase space structures that do not vanish in this limit will involve 
0{N) fermions at closely spaced energies. From the Young diagram perspective, 
0{N) rows with similar row length are necessary to give semiclassical structures. 

^°Well-localized Wigner distributions exist for coherent states which have a clean semiclassical 
limit. These distributions describe the circular motion in phase space that we would intuitively 
expect for a classical oscillator. All of this is in analogy with the fact that plane waves are completely 
delocalized in momentum whereas Gaussian wavepackets behave like classical particles. 
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The Fermi Sea: Consider placing the fermions in the lowest energy levels. 
This corresponds to the vacuum of the theory. There is no Young diagram describing 
the excitations because none of the fermions are excited (r„ = Vn in the notation 
introduced above). In the gravity description this should be empty AdSs x 5*^, which 
is obtain by the boundary condition M(0,r^) = 1 for < Ef and u(0,r^) = for 
> Ef. In the N ^ oo limit the Wigner distribution is 

oo oo 

2nhW::, = 2nhY,Wn{q,p) = 26"^^/^ J](-l)'^L„(4i7/?i) = 1, (126) 

n=0 n=0 

where H = {jp' + 0^)121 and we have used the summation formula 

Ln{x)z- = (1 - z)-' exp . (127) 

n>0 ^ 

We get the expected constant phase density for the vacuum, although the Fermi level 
has been pushed to infinity. Taking ^(0, r^) = Wsca. = 1 to be the boundary condition 
in the (xi, 0:2) plane does not lead to AdSs x S*^ as a solution because the asymptotics 
of the geometry are changed. For any finite A^ 

N-l 

27rhWse. = 27ihY,Wn{q,p) (128) 

n=0 

we find a function that has small oscillations around 1 for (g^ +j9^)/2 < hN and that 
approaches zero exponentially for larger E. A more detailed discussion of this can be 
found in [IH] ; here we will content ourselves by illustrating the point numerically. The 
scale for the decay and the oscillations is set by h. In the strict limit h ^ with hN 
fixed, 2nhWsea will be exactly one within the Fermi disk. Thus, in this limit we should 
map M(0,r^) = 27rhWsea{E) where = 2E. This is the familiar example of how the 
presence of many fermions with adjacent energies produces a smooth semiclassical 
distribution. We are also seeing that such a configuration in field theory leads to 
a gravitational description with a smooth semiclassical geometry. However, for any 
finite A^ and h the Wigner distribution contains data on h scale which maps onto the 
ip scale in gravity and is not refiected in the classical vacuum geometry. 

^^We thank Simon Ross for a discussfon concerning this point. 
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Figure 2: The sum of Wigner distribution functions of single fermions with excitation 
levels between n = 200 and n = 300 as a function of H/ h. 

The rectangular diagram: Another state in which many fermions have closely 
spaced energies is described by the rectangular diagram, e.g., 



(129) 



In such a state, each fermion has been excited by precisely the same amount 6 so that 
the operator creating the state has conformal dimension A = N6. In order for the 
energy 6 to be semiclassically visible we must take 5h finite as h —>■ 0. Since this is the 
same scaling as in the semiclassical limit, we see that A ^ A^^. This is satisfying 
because large classical objects in AdSs always have such large energies. Another way 
of saying this is that in the semiclassical limit we rescale both dimensions of a Young 
diagram by 1/A^ as A^ ^ oo and only those diagrams that survive this limit have a 
semiclassical description. The Wigner distribution corresponding to ()129|1 is 

NS NS 5-1 

27ihW,,,t = 27rhY,Wn{q,p) = 27r;i^ iy„(g,p) - 27rn^ (130) 

n=S n=0 n=0 

Following the discussion of the Fermi sea, in the strict ^ — > 0, Nh fixed limit, this will 
give a ring (for 6h < (p2 + q^)/2 < {N + l)6h) inside which the Wi gner distribution 
is equal to 1. Thus setting m(0;xi,X2) = W^ect{(l,p) in this limit reproduces the 
"black ring" boundary conditions described in However, for any finite A^ there 
are oscillations at the h scale both inside the inner radius of the ring and within 
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the ring itself (see Fig. Again it is clear, that the Wigner distribution cannot be 
directly mapped into gravity and a prescription is needed for removing the quantum 
oscillations. 



The triangular diagram: 

6 between any two fermions: 



Consider a state characterized by a fixed gap of energy 



(131) 



As discussed in Section 3.4, this is a typical state in an ensemble where the total 
number of D-branes is Nc and the conformal dimension is A = NNc/2. Triangular 
Young diagrams also provide the microscopic description of the superstar geometry 
[S], as will show in detail in Section In the N —>■ oo limit, the Wigner distribution 
of the triangular diagram is 



2nhJ2Wks{q,p) = 2e-'''/^J2^-l)''Lks{iH/h) (132) 



fc>0 



k=0 



with H = (g^+p^)/2, the energy in phase space. Next, summing ()127|) over the roots 
of unity {uj^,uj^, ■ ■ ■ ,uj^^~^^^^} with to = exp{2TTi/5) gives 



1 



fc>0 



j=0 



AH 00^ 



(133) 



Notice that the real part of uj^ /{I + uj^) is always 1/2. Therefore all terms have the 
same overall exp{2H/h) behavior. The terms with j oscillate, over a distance of 
AH ^0{h). Thus 



1 
~5 



27r/iiy^'^j^ngig = T + oscillations at scale AH ~ h 



(134) 



Even in the N oo limit, if h is finite, the oscillations persist, emphasizing that 
quantum oscillations in the Wigner distribution function cannot be erased simply by 
including many fermions. This is different from the familiar case of the Fermi sea 
(I126|) where the infinite sum erases quantum oscillations. Notice however that if we 
coarse-grain this density over a distance larger than the typical wavelength of the 
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Figure 3: The sum of Wigner distribution functions of single fermions with excitation 
levels between n = 200 and n = 300 separated by steps 5 = 2 as a function of H/ h. 



oscillations, all terms with j 7^ will average to zero, and what remains is only the 
term with j = 0. We find a coarse-grained average density of 1/6. 

Operators of any finite conformal dimension will of course correspond to diagrams 
for which the sum ()132|) is truncated: 



Numerically, one verifies that this Wigner function approaches a homogeneous ring, 
with superimposed oscillations, that stretches from H = Mq to Mq + M with density 
u = 1/6, so long as M is sufficiently large. Fig. El plots the same region as in Fig. |2 
but with 6 = 2. We see that the distribution in the classical region has its density 
reduced by precisely this factor 6. 

The examples above show how the large- limit serves to build up macroscopic 
distributions of fermions that can have a semiclassical interpretation. Nevertheless, 
at any the phase space distribution has important quantum mechanical features 
such as oscillations, and excursions below zero and above one. Since h in the field 
theory maps into ip in the bulk, we will argue in Section 4.3 that this structure refiects 
an underlying quantum mechanical structure in the spacetime at the ip scale. But 
half-BPS spacetime solutions of JH] are given in a supergravity theory that is effective 
in a semiclassical limit that coarse-grains over distances much larger than the Planck 
scale. Therefore, as a first step towards this semiclassical limit, we will study the 
structure of the Wigner distribution coarse-grained at the h scale and already find a 
structure that looks considerably more classical than the Wigner distribution, though 
it is still a complete quantum description of the system. 



M/S 




(135) 
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4.2.2 Husimi distribution 



Given the h scale oscillatory behavior of the Wigner distribution function, it is nat- 
ural to look for a second distribution function less sensitive to quantum mechanical 
structures, and more appropriate to the study of semiclassical physics. One way of 
achieving this is by convolving the Wigner distribution itself with a Gaussian kernel 
of variance h to smooth out quantum oscillations. For the effective single particle 
Wigner function W{p, q) ()12()j) we have 

Hu(g,p) = ^jdQdPZ{q,p-Q,P)W{Q,P), (136) 



Z{q,p;Q,P) = exp 



~{iQ-Qr+ip-pr) 



(137) 



The distribution function obtained in this way is called the Husimi distribution 
Hu„(g,p) jlH] and has several other useful meanings. For example, for particles 
in a harmonic potential, the Husimi distribution of a single particle state is a 
projection onto the coherent states \z): 

Hu^(g,p) = (138) 

{z\z) 



\z) = e^'^'-^'^lO) = e-l^l'/2 (139) 



where z = {p + iq)/V2h, \n) is the nth excited state of the oscillator, and {a,a^} 
are the standard ladder operators. For the eigenstates of the harmonic oscillator, we 
have the Husimi distribution 

Hu„(z) = K^|n)r = ^le-l^lVP"- (140) 

This is positive definite in contrast to the corresponding Wigner distribution ( p24j) 
and Fig. P). It has a single maximum achieved at |zp = n, which coincides with 
the classical orbit of the n-th eigenstate since {z]"^ = ijp' + q^)/{2h). (See Fig. HJ) 
The distribution function ()14n|l associated to single particle eigenstates is also ob- 
tained as the wavefunction in holomorphic quantization of the oscillator and as the 
wavefunction in the associated quantum hall problem j2SlEllllZ|- 

Finally, just as the Wigner distribution function computes expectation values of 
Weyl ordered operators, the Husimi distribution computes those of reverse normal 
ordered operators AM^q.p) in a state \ip) through 

(7/'|ijv(g,p)|V') = j dqdp Hu^(g,p)yl(g,p) . (141) 
This can be shown by using the identities 

(p-ig)"(V^|e(^'-'«)"^|0) = (7/^|a"e(P-'^)">) , (142) 
(p + ig)"(0|e(P+''?)"|V^) = (0|e(P+'«)"(at)"|V') , (143) 
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Figure 4: The Husimi distribution function of a single fermion with excitation level 
n = 50 as a function of H/ h. Compare with figure ^ 



from which it follows that 

J dpdqRn4q,p){p-tqr{p + tqr = {^\a"'{aY\^). (144) 

establishing the relation between the Husimi distribution and (reverse) normal or- 
dered expectation values. 

To study how peaked the maximum in the Husimi distribution is, we shall ap- 
proximate the full distribution by a Gaussian, and compute its width as a function 
of the excitation level n. Working with the variable C = (p^ + Q'^)/(2^), it is easy to 
show that for a single fermion the saddle point approximation gives 

Hu„(C) = ^ ^ e"^ C ~ ^ e-(^-")'/2n _ (145) 

The peak value is Hun(C = = n) = (2^7^- ^^^o Var(C — n) ~ and 
so increases linearly with energy; the standard deviation grows as ^/n. In terms 
of the physical energy variable E = (p^ + q^)/2 = h(, the standard deviation is 
(t(£') ~ ^/nh ~ y/Eh. In terms of the radial variable in phase space, r = a/p^ + = 
^/2E this translates into a standard deviation a(r) ~ ^/h. In this sense, the Husimi 
distribution gives a well-localized description of the harmonic oscillator eigenstates in 
terms of rings in phase space whose width shrinks as ^ — > 0. However, note that the 
spread in the variable E, a{E) ~ \/Eh is much greater than the separation between 
adjacent energy levels AE = h. In this sense, the Husimi distribution in phase space 
is still very far from being a ring localized in a thin band at the classical energy — 
the classical ring approximation is simply wrong for these states. While smearing the 
Wigner distribution with a Gaussian has removed the small scale oscillations, the h 
scale quantum uncertainties in the location and momentum of a single particle are 
still very much present. 
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We can also ask how the effective single particle Husimi distribution behaves 
when we assemble many fermions together as we must to obtain a semiclassical limit 
to compare with the half-BPS spacetime geometries. The analysis here is much easier 
than it was for the Wigner distribution because the answer is simply a superposition 
of the well-localized functions ()140|) . To assess how these superpositions behave let 
us study the N —>■ oo limit of the Husimi distribution for the filled Fermi sea and the 
triangular partition (jl31|) . Using the identity 

oo ^ (5—1 

^75;^^"")'"= 5 5:-"^^^ ' -=-"^'' (146) 

n=0 ^ k=0 

for (5 = 1, 2, . . ., it is easy to show that 
27rn Hu-.,,,, = 2nh J] Hu„,(C) = ] e-^^'--^^-^'^^^ ; C = (147) 

n=0 fc=0 

Thus in the N ^ 00 limit, the filled Fermi sea has a distribution 

27rnHusca = 27rnHu^ = 1 . (148) 
A triangular partition with row lengths increasing by 5 — 1 with 6 large leads to 



27r^ Hu^ia^ngle ~ ^ 



1 + 2e-^"''^'/''' cos(27rC/5) + ■■■], (149) 



where the additional terms are further suppressed by exponentials of (. It is easy 
to show that at the origin the Husimi distribution for a triangular partition always 
approaches 1, and decays exponentially to a stable value of 1/5 as C increases. The 
rate of decay is controlled by 6, but for any 6 of 0(1), the exponentially suppressed 
terms are irrelevant for the ( of 0{N) that are semiclassically relevant. At any finite 
and h, cutting off the sums at a finite limit produces exponential tails at the 
edges of the distribution that are determined by the variance in p45|) . However, the 
oscillations that were present in the Wigner distribution even after summing over 
many fermions have been smoothed out (see Figs. EI)- In all of these distributions, 
the relatively large variance displayed in (jl45|) is important for producing a smooth 
distribution from the sum over individual levels. 



4.2.3 Semiclassical Limit: the Grayscale distribution 

We will now systematically derive the semiclassical limit of the Wigner and Husimi 
distributions that should have a relation to classical half-BPS spacetimes. As we have 
described, the semiclassical limit requires taking ^ with hN held constant^^. 

^^For a recent discussion on the semiclassical limit of Wigner distribution function in the context 
of non- minimal strings see |48j . 



43 



Figure 5: The sum of Husimi distribution functions of single fermions with excitation 
levels between n = 200 and n = 300 as a function of H/h. Compare with figure [21 
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Figure 6: The sum of Husimi distribution functions for a triangular Young diagram; 
to = Nc/N = 300/100 = 3. The distribution plateaus at a grayscale of 1/(1 + cj) = 
1/4. 
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In such a limit the Fermi energy, namely the energy of the highest fermion in the 
vacuum, is kept fixed. Thus coordinates in phase space are effectively rescaled as 
q ^— q\/h and p i— >• p\/h. In the language of integer partitions, the states surviving 
the semiclassical large- limit are those for which rescaling the lengths of rows and 
columns of the Young diagrams by h leads to a finite limiting shape. In addition, 
quantum mechanical oscillations at the h scale should be averaged over; this has 
been partially achieved by the Gaussian convolution of the Wigner distribution that 
led to the Husimi distribution. To systematically derive a semiclassical distribution 
function on phase space, we should consider how semiclassical observables, namely 
operators that probe scales much larger than will respond in a given quantum 
state. One might worry that the ordering prescription for the operator will have 
a significant effect since the Wigner and Husimi distributions, computing the Weyl 
and normal ordered expectation values respectively, look so different from each other. 
However, in the /?. ^ limit the difference between ordering prescriptions vanishes: 
the semiclassical limit is universal. Hence we will simply analyze the semiclassical 
limit of the Husimi distribution, as it is somewhat easier to handle. 

We seek a new distribution function that is sufficient to describe the effective 
response of coarse-grained semiclassical observables in states that have a limit as 
^ — > 0. For the single Young diagram states of interest here this distribution on 
phase space, which we name the grayscale distribution, will be only a function of the 
energy E = {p^ + ?^)/2- To derive this, let AE be a coarse-graining scale such that 
/\E/h ^ oo in the h ^ Q limit. This simply means that AE must vanish slower than 
h in the semiclassical limit. A coarse-grained observable making measurements at an 
energy E responds to the integrated fermion distribution in phase space between E 
and E + AE. This response is given by 

pE+AE 

R{E,AE)=2nh dqdp'Ru{q,p), (150) 

Je 

where we have multiplied by 27rh to account for the overall normalization of the 
Husimi distribution. As discussed in the previous section, for a fermion of any fixed 
energy in the h ^ limit, the Husimi distribution becomes a thin ring. Therefore 
R{E, AE) counts the number of fermions with energies between E and E + AE. To 
obtain the effective density we have to divide by the area in phase space in this region: 

pE+AE 

Area= / dqdp = 2TiAE (151) 
Je 

where we used E = {jP' + q^)/2. Dividing ()150|) by p5H) gives the effective distribution 
that controls the response of observables coarse-grained at a scale AE: 

g{E) = 27rh \^^^^] • (152) 
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We can relate ()152|) to the Young diagrams describing specific half-BPS states by 
recalling that the number of fermions with energies between E and AE is simply 
given by the number of rows of the Young diagram with the appropriate number of 
boxes. In terms of the coordinate x introduced in Section 3 to index the rows of 
Young diagrams in the large- iV limit, 

^ Ax h 

Now consider any half-BPS state whose associated diagram is described by a limit 
curve y{x) as obtained in Section 3. Such a limit curve must exist for the state to 
have a well-defined semiclassical limit. Otherwise, as /i — > the structures in phase 
space become vanishingly small. Given such a curve, the energy of the xth fermion is 

E = h{x + y{x)) (154) 

Putting this together with p53j) gives 

9iE) = (155) 

where the right hand side is understood implicitly as a function of E computed by 
inverting (jl54j) to find a; as a function of E. This grayscale distribution captures the 
effective response of all coarse-grained semiclassical observables in the given quantum 
state. Remarkably, this is precisely the quantity ()11H) that we proposed on general 
grounds to determine the classical supergravity half-BPS solution associated to a 
Young diagram. Thus, explicitly, our proposal is 

u{0-r^)=giry2) = -^ (156) 
l + y' 

for all half-BPS states that are described by single Young diagrams in the semiclassical 
limit. To reiterate, we should first determine x as a function of r = ^/2E from ()154j) 
and then substitute that function into ()156|) to get the final form of m(0; r^). 



4.3 Quantum foam and effective singularity 

As we have discussed, half-BPS spacetime geometries are completely specified by the 
values taken by a certain function u{ri]Xi,X2) on the y = plane. Non-singular, 
topologically complex geometries result when m(0;xi,X2) assumes the values and 
1. Furthermore, an area quantization condition exists that requires that compact 
droplets in which u = 1 have their area quantized in integral multiples of 27r£p which 
was identified with h in the field theory. These facts, along with some explicit ex- 
amples, suggested an identification between the boundary values of the u function in 
the y = plane and the distribution of fermions in single-particle phase plane of the 
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dual fermionic picture of half-BPS states. Note, however, that the classical solutions 
of jlSI should be regarded as parametrizing a configuration space that still requires 
quantization. In particular in situations where the geometry displays Planck scale 
features which have high curvatures such a quantization is essential. An example 
is a single droplet of unit area 27r£p. Our analysis has used the dual field theory to 
find the exact quantum mechanical structures describing each half-BPS state. In this 
section we will argue that the exact structures imply an underlying "foam" at the 
Planck scale whose effective description in low-energy gravity is typically via various 
singular geometries. The singularity emerges because underlying quantum structures 
have been integrated out. 

Wigner and "quantum foam": First, suppose we map the Wigner distribution 
W associated to a half-BPS state directly to gravity as m = 2TThW. In general 
this function is bounded, but can be negative and certainly is not piecewise and 
1. Thus, a naive map of the Wigner distribution into classical gravity produces a 
singular geometry, and, following [321 EI] , the negative regions of W lead to tiny local 
closed timelike curves. However, both the period of oscillations in W and the size of 
the regions where the function is negative, are set in field theory by the h scale and 
so translate in gravity into the Planck scale ip. Thus the naive mapping from field 
theory to gravity produces singularities at a scale where classical gravity should not 
be trusted. A better interpretation is that the distribution is describing an underlying 
exact quantum structure that is relevant for the measurements made by certain kind 
of observables in quantum gravity that are related to Weyl-ordered operators in field 
theory. This underlying structure itself may not even have a good description as an 
actual geometry. A typical operator, encoded in a particular partition, has a detailed 
underlying Wigner distribution with complex oscillations. In the fermionic description 
the oscillations have two origins: (a) a regular oscillation associated to each fermion, 
and (b) randomized relative phases between contributions from different fermions, 
arising from the random spacings between fermion energies. This is describing an 
underlying "quantum foam" of extremely complex structures at the Planck scale that 
do not have literal geometric descriptions. 

Husimi and "classical foam": An alternative proposal would be to map the 
Husimi distribution Hu directly into gravity. This does not have quantum oscillations 
and is not negative, and builds up into coherent non-oscillating lumps when many 
fermions are piled up next to each other. However, in general it is not equal to or 
1, as we have demonstrated in the previous section. Thus, the effective geometry is 
singular. A typical state corresponding to a particular partition would be a complex 
superposition of the Husimi lumps in Fig. 01 spaced apart in a random way. If we 

Similar comments apply to the non-sing ular solutions of [7| El [TTl 1121 in the D1-D5 setting. 
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artificially approximate these Husimi rings as Is and Os in thin rings, this will give 
a non-singular effective geometry with complex topological structures (different S^s 
expanding) spaced randomly in the X1-X2 plane. We might call this a classical foam. 
However, because the typical structures will be at the h scale in the field theory, they 
will appear at the Planck scale in the geometry. Hence, the curvatures will be large, 
and we cannot really trust the classical picture. A better picture is that the non- 
singular geometries are a configuration space that should be quantized. The exact 
Husimi distribution then describes the underlying quantum mechanical object which 
is being approximated by classical geometry. 

We are in effect studying the quantum extension of the classical geometric con- 
siderations of and [H]. One might ask whether the Wigner or Husimi distribution 
is the correct extension, particularly given that they look so different. However, the 
distributions in fact give complete descriptions of the same quantum wavefunction. 
One is simply the natural description of measurements made by Weyl ordered observ- 
ables, while the other describes reverse normal ordering. Indeed, each distribution 
can be transformed into the other. The Husimi distribution is the convolution of a 
Gaussian with the Wigner distribution 



where we have Fourier transformed the phase space variables q and p to turn the con- 
volution ()157p into a product. Using the definition of correlation functions ()116p and 
1)1411) with these relations, it is easy to show that reverse normal ordered operators 
that are non-local in phase space can recover that same information from the Husimi 
distribution that is encoded about Weyl-ordered operators encoded in the Wigner dis- 
tribution. The lesson about gravity is that the underlying "foam" looks different to 
different observables. Local Weyl-observables see a more quantum mechanical struc- 
ture, while local reverse-normal-observables see a more classical structure. This sort of 
phenomenon has been seen before in string theory — D-branes, strings and gravitons 
can propagate in rather different effective geometries within the same spacetime. 

Effective singularities: Of course a description in classical gravity presupposes 
that the semiclassical limit has been taken and that we are looking at observables 
that probe lengths much larger than the Planck length. In the previous section 
we computed the phase space distribution ()156|) appropriate for such semiclassical 
observables for states that exist in the h ^ limit. We argued that this effective 
grayscale distribution is universal because the difference between Weyl and normal 




(157) 



and this can be inverted as 
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ordering prescriptions vanishes in the h ^ hmit. One potential subtlety is that 
for operators of high dimension, commutators of many creation and annihilation 
operators can give rise to c-numbers of order N, in which case there might be a finite 
correction even in the classical limit. We will ignore this interesting possibility for 
the present. We found that 



It is clear that states that do not correspond to rectangular, or step, diagrams {y' ^ 
0, oo) give boundary conditions in gravity that will lead to singular geometries. The 
full geometry is constructed in terms of a function 



This equation is obtained from ()100|) for rotationally invariant functions z{0; r^) by 
performing the angular integral and where the integration variable t is related to 
r', the original radial variable in phase space, hj t = r'^ + 1]"^ — r^. This gives the 
universal low-energy spacetime description of quantum states described by Young 
diagrams that lie close to the limit curve y{x). In this way, a singularity in spacetime 
arises because underlying quantum structure at the Planck scale has been integrated 
out. 

What information is lost in this semiclassical limit? There are three obvious ways 
in which the description in classical geometry loses quantum detail: 

1. The effective description of all states that lie close to a particular limit curve 
is the same. Thus the detailed difference between these states has simply been 
lost. For example, since almost all half-BPS states with a fixed total charge A lie 
close to the limit curve PHj) . almost all such states will have the same description 
in classical gravity. The differences between these states are encoded in gravity 
in Planck scale structures. 

2. Given a particular microstate, the grayscale distribution appropriate to semi- 
classical observables ()156p erases quantum details of the exact distribution on 
phase space. These details are encoded in gravity in Planck scale structures that 
do not appear in an effective long-wavelength description in classical gravity. 
The relevant Planck structures do not have a description in terms of smooth 
geometries. 

3. The half-BPS quantum states are encoded in terms of a 2A^-dimensional A^- 
fermion phase space, but only the effective two-dimensional one-particle phase 
space makes an appearance in the geometric description. Some quantum ob- 
servables, such as massive string states, could be sensitive to details of the full 
A^-particle phase space which has been erased. 



1 



1 



(159) 





(160) 
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We will return to (3) in Section 5. Below we will discuss the explicit construction 
of the effective geometries describing the generic half-BPS states of fixed charge and 
number of D-branes (which includes the well-known superstar jHI), and the generic 
half-BPS states of fixed charge and unrestricted number of D-branes (a new space 
that we will call the hyperstar). 

The Hyperstar: Using the limit shape ()43|) to compute the slope y'{x), we can 
derive the grayscale by substituting into (jlllj) : 

z{0; r^) = ^-t g^(^) ^ u{0; r^) = t g^(^) = 1 - g^"^, (161) 

where t = = 1 — g^. In order to determine the full geometry, we need to 

relate the continuous variable x with the radial variable r^. This can be done by using 
the identity y{x) + x = r'^/{2h), which we rewrite as 



r 



2 



log (1 - g^-^) log (g^ + t qP) 



P= — = C{I3,N)+ "\ '- + X ^x=^^ ^. (162) 

in logg logg 

The grayscale distribution in phase space is then 

which is nothing but a standard Fermi-Dirac distribution with chemical potential 

= q-N -I (164) 



also considered in The full type IIB configuration is therefore characterized by 
the function: 



O - y y^^_^^ dt -^—-^—^^^-^-^ [t2 + (2r77)2]3/2 • 

Just as the limit shape PHj) has different behavior in different regimes, the same 
fact applies to the distribution function ()163p . Consider the bottom of the Young 
diagram, that is, (x/N) <C 1. Since the limit shape is linear in this region, we 
know the dominant contribution to the distribution function should be a constant. 
We shall also compute the first correction to this. In this regime, the distribution 
function written in phase space variables is 

u{0- r') = l-e~--a^{l- e"") ^ . (166) 

There is an upper bound on the radial coordinate due to the condition a; -C 
under which these expressions were derived. This bound satisfies 
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Thus, as expected, this is a regime in which we are describing the phase space dis- 
tribution at radial distances smaller than the scale associated with the Fermi surface 



of the vacuum [s/hN). We can replace with t and integrate this to determine the 
LLM function z{ri; r^) in this regime. 

Near the top of the Young diagram, that is, at a; ~ A^, the grayscale distribution 
is approximated by 



a 



X 



X 



N 



(168) 



In this regime, the limit shape is a logarithmic curve, from which we can extract the 
relation between the radial coordinate in phase space and the row coordinate x: 



X 

N 



1--W 

a 



-ar^/2Nh 



e° 1 



(169) 



where the right hand side is written in terms of Lambert's VT-function, which is ob- 
tained as the inverse function of the equation W = f{W). The droplet distribution 
is 

-ar'^/2Nh 



W 



1 



(170) 



This regime probes scales in phase space much more energetic than the Fermi surface. 

We should note an important consistency check to our analysis. We have asserted 
that the conformal dimension is 



A 



— - — M(0:r2) - - 
2'Kh 2h ^ ' ' 2 



(171) 



which reproduces the expected results for semiclassical geometries which are non- 
singular. This also applies to singular spacetimes like the hyperstar. As the top row 
of the limit curve is infinitely long, we must integrate the excitation level p of the 
fermions out to oo. We have 



A 



log(l + tq 



.p~N\ 



P- 



\ogq 

The second term gives 



+ 



Lis 



J 



(logg)2 

log(g^ + 1) 



\og{q + tq 
\ogq 



N 



logg 



(172) 



(173) 



As this term should fix the number of fermions (we are subtracting out the ground 
state energy corresponding to empty AdSs x S*^), we see that t = 1 — q^ . Comparing 
to (jl21), we have indeed correctly reproduced C{P,N). Evaluating the first term as 
well, we now have 



Li2(l - 
(logg)^ 



-N\ 



2 



Li2(l - q 
{log qy 



(174) 
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e 



0.-_q-iN^Na)^-l^^_^^^S^ ^. (177) 



where we have used a dilogarithm identity Li2{z) + Li2(^3Y) I ^'-'^(l ~ -2)^ = 0. This 
is the expected result. 

The ensemble containing the Superstar: For a generic temperature, the limit 
shape describing a typical state in an ensemble in which the number of columns (D- 
branes) was fixed (iVc) is given by equation (jH^ . We can compute the slope y'{x) of 
this curve at a given point, and from it, using our proposal (jlllj) . derive the phase 
space distribution 

n - (1 - C)g^(^)-^- (within droplet), 
^ ' ' \ I (outside droplet). ^ ' 

Using the identity y{x) + x = r'^/2h and the limit shape (j84j) . we then invert x as a 
function of r^: 

I +,^.,ry,n ' (176) 
where e^^ again plays the role of a chemical potential, 

-N 1 
l~qNc 

The corresponding full type IIB configuration is determined by the integral expression: 

^^'^ 2 y^2_,2 «^ i + e/3/.e-/3(t+r2-^2)/2;j [^2 + (2rr7)2]3/2 • 

This assumes the same functional form as the corresponding equation for the hyper- 
star ()165|) : the differences are in the definition of the chemical potential and in the 
upper bound of the integration. The upper bound is necessary because the droplet is 
compact. There is no infinite tail to the Young diagram as in the hyperstar case. 

As discussed in Section 3.4, in the infinite temperature limit the typical state is 
described by a triangular Young diagram. In the large- limit this is described by 
a straight line limit curve (jHOJ. The grayscale is a constant. Using (jllip . the phase 
space distribution zsiO-jv"^) will also be constant: 

• ' \i if rV2fi > /V + /Vc. ^ ' 

Since, within the droplet region, this number is different from ±1/2, the spacetime is 
again singular. 



^'*Since the distribution u(0,r^) has noncompact support, there is some question as to whether it 
describes an asymptotically AdSs x geometry. Of course u falls off exponentially at large r^, but 
AdS geometries also expand rapidly at large distances. 
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In the following we will show that the geometry derived from ()179p is precisely 
that of the superstar jS]. This will be both a check of our formalism, since we have 
already physically motivated that our ensemble analysis is the appropriate one for 
describing the microstates of the superstar, and also a confirmation of the emergence 
of a singular spacetime as an effective description of microstates that differ from each 
other by Planck scale structures. 

Let us first remember that the energy (conformal dimension) and total number of 
D-branes of the superstar configuration j3] are given by 

A = w— , Nc = u:N , (180) 

where uj = qi/L"^ is the parameter describing the charge of the system. Notice that 
the slope of the limit shape (jnn|) is precisely equal to u, i.e. y{x) = ux. It is now 
clear that the area under the limit shape equals the conformal dimension (A) given 
above. Noting that 

MO;r^) = i-..(0;r^) = -^;-^, (181) 

in the region of the phase space plane between r'^/2h = and = + Nc, we 

may also readily verify from (jlTlj) that the distribution function also reproduces A. 
The quantized number of giant gravitons Nc is encoded in the radius of the gray 
droplet, or as the integral of (1 — us) over the droplet. 

Let us rewrite the ten-dimensional metric corresponding to the half-BPS superstar 
configuration in type I IB: 

ds^ = -^fdt' + ^dr'' + ^r^ds% + ^L''del + ^smH^dsl, 



+ — cos^ 01 \L d4>i + {H^^ - 1) dtf , (182) 

where Hi = 1 + gi/r^, / = 1 + Hi/ L"^ and 7 = 1 + gi sin^ Oi/r"^. If we compare the 
physical size of the two three-spheres that appear in the superstar metric with their 
parametrization in ? we obtain the conditions: 

r^e^ = v^r^ , 

r L"^ sin^6'i 
r/e-^ = 183 

Solving the system 



= r2 sin^ Qi , e« = , (184) 

L sm 61 

and using the fact that the grayscale distribution is given hy z = (l/2)tanhG' we 
obtain: 

1 r^7 - sin^ 61 



2 r'^'-^ + sin^ Oi 



;i85) 
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Since it is the value 2(77 = 0) that is related to the semi classical distribution function, 
we must analyze the behavior of G at rj = 0. We observe that there are two different 
coordinate regimes for which this applies: 



(i) When sin^i = 0, z{r] = 0) = 1/2. This is consistent with the fact that whenever 
sin^i vanishes, the distribution of D-branes vanishes. There are no fermion 
excitations in this locus. Since we follow the conventions in [T3], 1/2 is the 
right boundary value to describe absence of excitations in the fermionic picture. 

(ii) When r = 0, the density of giants is non-vanishing. We get 



It is reassuring to check that the expression ()186p identically matches what is derived 
from purely field theoretic and statistical mechanical considerations ()179|) . 

5 Correlation functions 

In Section 2.2 we argued that correlation functions computed in typical very heavy 
states are to great accuracy independent of both the operators that appear in the 
correlation function as well as the state that is being probed. Correlation functions 
therefore only depend on charges and other quantum numbers that clearly distinguish 
states and operators at the macroscopic level. Deviations from universality are ex- 
pected to be of the order of exp(— 5*), with 5* the entropy of the ensemble from which 
we pick the states that we probe. In this section we will try to make these claims 
somewhat more precise for the case of half-BPS states. 

In general, correlation functions computed in a half-BPS state will involve the full 
Yang-Mills theory because the intermediate states in an arbitrary correlation function 
can explore the full theory. However, there are special classes of correlation functions 
which can be computed entirely within the half-BPS matrix models. The authors of [31 
121] describe such correlators and group theory techniques for computing them. Here 
we will take a different perspective and simply use the Wigner (or Husimi) distribution 
function to compute correlation functions as phase space integrals. While we will 
will not compute the general correlation function of the theory, it is nevertheless 
interesting to ask whether these matrix model correlators are sensitive to the difference 
between typical states. In general we expect our results to be renormalized in the 
full theory, but the matrix model results give useful insights anyway. It would be 
interesting to work out the representation of the operators with exact matrix model 
correlators in the language of functions on phase space. 
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5.1 Phase space approach 

The natural operators in the half-BPS matrix model are products of tr(X^). Such 
operators can be rewritten in terms of tr(A^) and tY{{A^)'^), the matrix creation and 
annihilation operators. By definition, all correlation functions computed within the 
half-BPS matrix model can be equivalently calculated through phase space integrals 
of the corresponding phase space distribution functions. Generic multiple trace opera- 
tors involve 2A^-dimensional phase space integrals, and so they require the knowledge 
of the 2A^-dimensional phase space distribution. On the other hand, the semiclassical 
gravitational description of any typical state in the ensemble only depends on the 
2-dimensional single-particle phase space distribution function, raising the question 
of what data about the full quantum mechanical system survives the semiclassical 
limit. 

The special class of single-trace operators has expectation values that can be com- 
puted by just knowing the single-particle phase space distribution. Such operators, 
when expressed in terms of eigenvalues, become a sum ^) "with Aj the 

jth eigenvalue of X. This is a sum of one-fermion operators, and thus its expecta- 
tion value can be completely computed using the single-particle Husimi or Wigner 
distribution u{q,p). We will focus our attention on such operators here. 

We are interested in computing the expectation values of single trace operators 
in states appearing in certain ensembles, as well in measuring the variance over the 
ensemble of these responses. The variance to mean ratio will give a measure of the uni- 
versality of correlation functions computed in typical states. Let us denote by Ux{x) 
the Wigner or Husimi distribution of the partition A (state), where x = {q,p) stands 
for the phase space point. The ensemble average of the distribution function u^{x) 
in an ensemble £ with weight w{X) will be given by Uw{x) = X^^ef ^(-^)^^(^)- 
definition, the expectation values of observables F and their corresponding ensemble 
averages are 



respectively. On the other hand, the variance in this observable over the ensemble is 



Var^(F) = / dxdyF{x) F{y) 5^«;(A) [ux{x) - u^{x)] [ux{y) - u^iy)] . (190) 



Xe£ 

As expected, this is an expression that depends on the fluctuations of ux{x) around 
its mean ensemble value Uw{x) at different points x. 




(187) 



(188) 




(189) 



Using the definitions for Fx and {Fx)u,, this variance can be written as 
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For the subset of states considered in this paper, all phase space distributions are 
rotationally invariant, and so we write them as u{r), r being the radial variable in 
phase space. Furthermore, from our canonical ensemble analysis, the fluctuations in 
the number of Young diagram columns of different lengths are uncorrelated. It is easily 
shown that this translates in phase space into the independence of the fluctuations 
of ux{r) around its mean ensemble value M«,(r) at different points r. In other words, 

^ w{\) [uxir) - u^{r)] [ux{r') - u^{r')] = 5(r - r') ^ w{\) [ux{r) - u^{r')]'^ . 

(191) 

Using this, we can rewrite p90|) as 

Var^(F) = J dxF\r) [{uD^ - {ux)l] , (192) 

where we also assumed for simplicity that the operator F itself was rotationally 
invariant. Therefore we learn that the variance in observables over the ensemble £ is 
controlled by the variance of the phase space distribution over the same ensemble. 

Let us first consider very delocalized observables spread over the entire compact 
region where the distribution ux has support. (Of course ux may have exponential 
tails, but these will not make a large contribution, so we can ignore them and speak 
about the region of support as being compact.) Then for any reasonable observable 
F can be bounded by some Fq within the region of support and we can bound the 
mean and variance of the observable by 

F^ < Fo j dxu^{x) = FoN, (193) 

Var^(F) < F^ J dx {uD^ < F^ j dxu^x) = F^ N . (194) 

For sufficiently delocalized observables this leads to the estimate 

Var^(F) ^ ]_ 
Fl ^ N 



(195) 



so that deviations from universality behave like . For most ensembles we consider 
the entropy scales as ~ N, so this is in agreement with our general considerations 
in Section 2.2. 

Next consider a highly localized observable F{q,p) ~ Fq 5(A/gM-~p^ — tq). In fact, 
such a probe localized below the h scale is unrealistic, but it is instructive to study 
it anyway. In this case 

~ Fo{ux{ro))^, (196) 
Var^(F) ~ F^ [{ul{ro))^-{ux{ro))l]=F^VaTUux{ro)). (197) 
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In the semiclassical limit, both the Wigner and the Husimi distributions approach the 
grayscale distribution of Section 4.2.3 and so in this limit we take ■u^(r^) = 1/(1 + y') 
in terms of y{x), the Young diagram limit shape. Using standard error propagation, 
it easy to relate the variance in u for small r^, which is the region containing most of 
the structure of half-BPS states, to the variance in the y'{x). Doing this one shows 
using that 

Var^F) ^y^^^^y.^^2j^^^ (igg) 



Finally, by definition 



i?2 

w 



, In-(x+Ax) '^i In-x '^^ ^« ,^ „„s 

y = — ^ '—^ ~ cn-x (199) 

where cn-x is the number of columns of length — x in the Young diagram. From 
this, 

^-^(N-x) 

VaTu,{y'{x)) = VaT^{cN-x) = _ ^.^(^N-x)y (^00) 

Putting this together with ()198j) and recalling that f3 is 0{1/N) for the ensembles of 
interest we see that for small x (near the origin in phase space) where ~ 1, the 
delta function observable has 0(1) variations over the ensemble. For x of 0{N) (where 
the error propagation formula we used is not strictly valid) the enhancement of the 
variance by the denominator in ()200|1 competes with the exponential decay in Uw In 
any case, the main point is that extremely localized observables show larger variation 
over the ensemble and thus can be used to identify details of the state. However, the 
S function observable considered here will be extraordinarily difficult to construct as 
a probe in the dual gravity theory. Recall that the h scale in the field theory maps 
into the Planck scale in gravity, so we are describing the results of measurements 
by a super-Planckian probe. Such things will be extremely difficult to construct in 
gravity, because gravitons, strings, D-branes and black holes, which are the natural 
probes in string theory, all have sizes much bigger than the Planck length. Indeed, 
in general adding energy to a probe in string theory tends to increase its size. In 
fact, to construct probes that measure super-Planckian scales, non-local objects may 
be needed. To see this, observe that in the dual field theory that we are discussing 
here, the natural probes are polynomials in X with traces distributed to impose gauge 
invariance. These objects map into the usual gravitational probes. Local probes in the 
field theory simply cannot be localized below the h scale according to the Heisenberg 



Strictly speaking, in the semiclassical limit we should not be using observables that are localized 
at scales smaller than h. But we will do this anyway to get a sense of the effects of localizing a probe 
in phase space. 

^^In the small r^, the fluctuations are small, but there are enormously many configurations. At 
large the fluctuations are larger and the error propagation formula is not really valid, but there 
are also many fewer configurations to consider. 
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uncertainty principle. Thus, one is led to consider non-local combinations of objects 
involving arbitrarily large frequencies. 

A particularly efficient set of observables {Fn{q,p)} would be functions of phase 
space that are the single fermion distribution functions. This is because the full phase 
space distribution is simply a linear sum of such constituent terms. Thus integrating 
the Fn{q,p) against the distribution function for an individual diagram A acts as a 
filter for testing whether the associated state contains the excitation level n. In this 
way it would appear that there is an 0{N) set of observables that can efficiently 
"detect" the half-BPS state associated to a given Young diagram. This phenomenon 
is occurring because the half-BPS states are an integrable sector of the full Yang- Mills 
theory. This leads to the free nature of the system and corresponding simplicity of 
identifying states. However, note that the filtering observables {Fn{q,p)} will not be 
simple to construct in terms of the multi-trace observables that are naturally related 
to probes in gravity. We should expect that they are related to non-local observables 
of the Yang- Mills theory like other integrable charges of this system jlH] ■ Nevertheless, 
the basic idea here that certain probes extract patterns from the underlying quantum 
state is in analogy to the discussion of Schwarzschild black holes in Section 2. 

Another approach to computing correlation functions using the fermion phase 
space is given in Appendix A. 

Exact correlation functions in phase space: Since we know the exact single 
particle quantum mechanical Wigner distribution function for a given Young dia- 
gram, it is possible to compute the exact matrix model correlators for single particle 
operators. To be definite, let us consider some observable whose phase space de- 
scription is given in terms of p*g*, i.e. F{p,q;t, s) = p^ q'^ . We wish to compute 
the expectation of the Weyl ordered operator corresponding to g; t, s) in a state 
corresponding to a partition A. From Section 4, the Wigner distribution of the state 
is W\ = ^e~(P^+^^) J2mex ("1)"" -^m[2(p^ + g^)]. Working with polar coordinates in 
phase space, the expectation value is then given by 

poo /'2lT 

Fx{t,s)= / rdr / rf0 r*+" sin* cos' • (201) 
Jo Jo 

Due to rotational invariance of the state, the above integral is only non-vanishing 
when both (t, s) are even {t,s) = {2i,2n). Using this fact, and some identities, the 
result is as 

m 

(202) 

A similar exact expression can be given for normal ordered correlators using the 
Husimi distribution. From this point of view, the question of identifying the partition 
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A from the set of observables F{p, q; t, s) amounts to asking how much detail in ()202|) 
is needed to tell apart correlations in two different typical states. This might be 
a convenient formulation of the problem. It would also be useful to translate the 
functions F{p, q; t, s) back into the language of gauge invariant operators of the Yang- 
Mills theory in order to identify precisely what observables these are in string theory. 

Multi-particle phase space: If F is not a sum of one-fermion operators, for 
example if we are interested in computing the correlation function of tr(M^)tr(M') 
in the matrix model, the result can no longer be written as a simple integral over the 
single particle phase space. For ensembles like (jl87|) . one can show that expectation 
values of operators that can be written as sums of t-fermion operators reduce in 
general to t-fold integrals over phase space, with rather complicated kernels that are 
functions of the t-points. This is quite natural from the point of view of the AdS / CFT 
correspondence, where t-point functions involve multiple integrals over the bulk as 
well. In the present case the bulk is in general singular, but we nevertheless find 
concrete and precise expressions for the multipoint correlators as multiple integrals 
over phase space. Thus, in this half-BPS context, we have in principle a precise 
proposal for the boundary conditions that need to be imposed at the singularities. It 
would be interesting to explore this further. 

5.2 Second quantized approach 

We can try to improve upon the above analysis in a second quantized formalism. 
Consider a free Fermi field with a Hilbert space with modes ipm, m E Zi + ^, and with 
commutation relations {4'm,'^n} = ^m+n- We will restrict attention to the A^-fermion 
sector of this Hilbert space. The action of the matrix creation and annihilation 
operators is given by the following fermion bilinears 

n>0 

trmn - E (203) 

n>0 

and one can in a similar way work out more complicated single trace operators that 
consist of a product of terms and {A^y. 

To illustrate the computation of correlation functions, consider 

O = tT{AP)tT{{A^y). (204) 

For quantities like ()187|1 . we only need the piece of O that acts diagonally on the basis 
of states consisting of partitions A. If we insert ()203|) in O, the only contributions to 
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(jl87|) are therefore from 

{n + p)\ 



^0 = E -Z^^-hn^^,^p-,n^-h-P-n^h^n- (205) 



n>0 



Such an operator gets contributions only from fermionic states with both the n-th and 
(n + p)— th harmonic oscillator levels are occupied. When computing the ensemble 
average ()188p . it is natural to expect to get the same contributions as above, but 
weighted with the probability that such excitations occur in the ensemble. Since in 
this paper we have only considered thermal fermionic distributions, these probabilities 
should be given by the standard Fermi-Dirac distribution:^^ 

I + q"- 

Thus the ensemble average (jl88|) of O will be equal to 

= V ^-^^-^v{n)v{n + p). (207) 

n>0 

For small values of p this can be approximated as follows 

n>0 

fOO 

dnnFv{nf. (208) 







More generally, we can consider integrals J dn n"v{nY . If we expand the denominator 
in v{n) and do the integral, this reduces to 



a+l 



[ dnn-v{nf = r(« + l)(i^)" + / ) 

- r(a + i)f-:±y''V-i^g-^('3+^)(i + r)/^-" 



-1 o-A'C«-l) 



In particular, turning back to ()208|) we get 



_1 \ 



{O)^ ~ -V{p + 1) ( 1^ ) g"''Li,(g-^). (210) 



^^This is consistent with the fact that, in the semiclassical hmit, the phase space distribution 
obtained through our analysis of the hmit shape Young diagram is precisely given by the Fermi- 
Dirac distribution. 



60 



There are a couple of interesting features in this expression. First, its scahng with N 
is controlled by (— l/logg)"^"^. Second, as we claimed in Section 2.2, the ensemble 
average does only depend on the parameters specifying the ensemble {/i, /?}. These 
are in one-to-one correspondence with the global charges {A^, A} characterizing the 
set of microstates under study. Thus we see that once we average over the ensemble, 
which is the right way to infer the properties of correlation functions in typical states, 
the response is mainly determined by their global charges. 

Let us compute the variance in the expectation value of O to estimate the univer- 
sality of the above result as we consider different typical states. The only contributions 
to ()189|) from O are coming from the diagonal part of O which was explicitly written 
down in ()205|) and so it equals 

Var,(0) = {{Oo)') - {Oof. (211) 

Just as the operator Oq gets contributions if the harmonic oscillator eigenstates n and 
n + p are occupied, (Co)^ gets contributions if eigenstates n, n + p, m and m + p are 
occupied for some m, n. Generically these are four different eigenstates, but when \m— 
n\ = p there are only three and if |m — n| =0 only two. These special configurations 
of m, n are the only ones that contribute to 1)2111) in thermal ensembles. All other 
(generic) configurations of n, m are uncorrelated and cancel against each other in the 
two terms in 1)2111) . As in ()207j) . when we compute the ensemble variance, we weight 
these contributions with the probabilities for such harmonic potential excitations to 
occur in the ensemble. This way, the ensemble variance is 

^ 2w) ! 

Var^(O) = 2^^-^ -^v{n)v(n + p)v(n + 2p)(l - v(n + p)) 



n>0 



n>0 

In view of (12091) this behaves as 



+ [v{n)v{n^p)-v{nfv{n^pf). (212) 



~ pgi^ (213) 

and in particular 

Since logg behaves as 1/-\/A=1/N, we find the same behavior in the variance of the 
correlation function as we found in the previous example. This again agrees with 
the general expectation from Section 2.2 that the variance should behave as 1/ logS", 
with 5* the entropy. 

When does the variance become large? Superficially, the prefactor in ()214|) be- 
comes of order one for p ~ log A^. This is not quite what we would expect, we would 
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expect a breakdown for p ~ A^^ for some 7 rather than for p ~ log A^. This may well 
be due to the various approximations that we made in these estimates and we leave 
a more detailed study of this to future work. Amongst other things, general correla- 
tions cannot simply be computed in the matrix model, an integrable subsector of the 
Yang- Mills theory, as we are doing here. We can explore the point where the variance 
becomes of order one also emerge when we consider multi-trace operators that will 
probe k different harmonic oscillator levels at the same time. Then in computing the 
variance as we did above, there are roughly k"^ possibilities for special arrangements 
of the eigenvalues probed by the two operators where at least two of them coincide. 
Thus we obtain an enhancement of the variance by an extra factor of fc^, suggesting a 
relative variance of order one when k ~ N^/'^ . More analysis is needed to define and 
analyze the universality of correlation functions. 

6 Discussion 

In this paper we have shown that very heavy pure states of quantum gravity can have 
an underlying quantum description as a "foam" with complex structure at the Planck 
scale, the details of which are invisible to almost all probes. The "foam" might be 
more quantum mechanical or more classical depending upon the nature of the probe. 
Whether one has a useful geometric description of the underlying microstate de- 
pends on the observable that one is measuring. Specific classical geometries give the 
universal low-energy, long-wavelength descriptions of large classes of the underlying 
quantum states. We showed how these classical descriptions are derived and how 
smooth underlying quantum dynamics can give rise to singular effective semiclassical 
geometries. As an example, we explicitly constructed the singular spacetimes giving 
the universal low-energy description of almost all half-BPS states with a fixed mass. 
The striking resemblance between the reasoning for the half-BPS states and the quali- 
tative analysis of Schwarzschild black holes in Section 2 suggests that our exact results 
in the super symmetric scenario give a generally correct qualitative understanding of 
the non-supersymmetric case. 

In the situations we have described, it is explicit that the low-energy classical 
description erases many crucial quantum mechanical details that are present at the 
Planck scale. Information loss is thus implicit in the classical description of spacetime. 
But when the detailed quantum mechanical wavefunction is incorporated there is 
no room for a breakdown of unitarity. It has been suggested in |5U[ 1^ that the 
information loss paradox for black holes is solved because the wavefunction for a black 
hole spacetime also has support on geometries that do not have a horizon. However, it 
was also argued in [521 ESI ^^^^ summing over geometries in this way is insufficient to 
restore information apparently lost in black holes. We are proposing that the quantum 
mechanical wavefunction is indeed important, but sums over particular smooth and 
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singular geometries using the Euclidean approach to quantum gravity will not suffice 
to restore unitarity. The black hole in our picture is simply the universal effective 
low-energy description of a complex underlying state. The loss of information lies 
precisely in the universality of the semiclassical description.^® 

To see an actual horizon emerge in a low-energy effective description, one will have 
to consider states will less supersymmetry. It may also be that interaction effects at 
finite coupling will become important. Any star will collapse to a black hole if the 
Newton coupling is large enough, indicating the appearance of horizons does depend 
on the coupling. This is despite the fact that Qs drops out of the expression for the 
conformal dimension of a black hole operator. It will be important to understand 
this issue and whether there is any phase transition as a function of the coupling 
involved. Nevertheless, it seems that the heavy states that we are describing know 
about gravity and don't just behave like thermal gases. Probes of a standard thermal 
gas in a field theory that have energies above the thermal scale pass through the 
gas completely untouched. In our case, a probe of very high conformal dimension 
interacts very strongly with any very heavy state simply because there are more ways 
for the many fields in the state and probe operators to interact. We suggest that this 
is an essentially gravitational phenomenon — heavy objects interact more strongly 
because, in gravity, the mass is the coupling constant. 

Our results raise many questions including the following: 

1. The limit shapes we have discussed give a "mean field" description of the mi- 
croscopic physics. The limit shape translates to a semiclassical geometry on the 
gravity side and gives universal correlation functions in the field theory because 
small-scale fluctuations have been integrated out. In a sense, the limit shape 
is a "master field". Fluctuations about the limit shape should be thought of 
in terms of continuum field theory. It would be interesting to understand this 
structure more explicitly. 

2. Can we assess correlation functions in greater detail to understand what infor- 
mation about the microstates of gravity can be extracted by different kinds of 
probes in both the half-BPS sector, and more generally without supersymme- 
try? We expect that information about global charges is easy to determine and 
that detailed information about the nature of microstates will require carefully 
designed probes capable of making exceedingly fine measurements. The ele- 
mentary analysis of the Schwarzschild black hole suggests that it is far easier to 
decide what a black hole is not rather than what it is. Distinguishing the state 

^^It seems to us that recent attempts to "see behind a horizon" using AdS/CFT [SJ will also not 
be able to identify the state underlying a black hole precisely because these efforts never account for 
the details of particular microstates. It would be interesting to understand whether the structures 
in the complex time plane appearing in these papers arise in some effective way from the underlying 
"foams" described in this paper. 
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would at least require determining all of the multipole moments. How may we 
do this efficiently? 

3. We have discussed a basis of half-BPS states related to single Young diagrams 
preserving U{1) symmetry. Generic half-BPS states, which are sums of these 
diagrams break this symmetry. A particularly convenient method for describing 
the semiclassical limit of such states is in terms of coherent states. This gives 
rise to a semiclassical picture of topological complex Planck scale droplets whose 
coherent sum gives rise to smooth geometries. What can this teach us about 
the underlying quantum structure of gravity and the emergence of classical 
spacetimes? 

4. While the N fermions that appear in the description of half-BPS states have 
a 2A^-dimensional phase space, only the effective two-dimensional one-particle 
phase space makes an appearance in classical gravity. What is the role of the 
full A^-particle phase space and what does the full quantum system in the field 
theory teach us about boundary conditions at the singularity in the effective 
classical spacetime? 

5. How do the lessons from the analysis of typical half-BPS states in AdSs x 
translate to other settings? Can a similar picture of underlying "foam" and 
effective singular geometries be developed, for example, for the D1-D5 string 
that underlies most extremal black holes in string theory? Some efforts in this 
direction, using the techniques of Mathur and collaborators 0, will be reported 
in [SSj. 

6. Can stringy corrections "cloak" the singular effective geometries describing our 
generic half-BPS states with a horizon, in a manner similar to ISEI? Does the 
area of the horizon match the entropy of the associated microstates? 

7. What is the relation between our analysis of typical half-BPS states and the 
description of quantum foam in the topological string jSZI? In particular, there 
are strong similarities in the appearance of limit shapes of Young diagrams 
and topologically complex geometries containing random insertions of blown- 
up spheres in an otherwise smooth and locally flat manifold. 

8. In Section 2 we gave arguments justifying our expectations that the conformal 
dimension of heavy operators (A ~ A^^) would not renormalize strongly. It 
would be interesting to analyze this point precisely. Along the same lines of 
extrapolating from weak to strong coupling, it would be interesting to under- 
stand better the role of (possible) phase transitions in the appearance of bulk 
horizons at strong coupling in the dual field theory formulation. 
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9. In principle, there is no upper bound to the dimension of operators in the Yang- 
Mills gauge theory. What do typical states of conformal dimension greater than 
A^^ look like in gravity? 

10. In the phase space droplets giving rise to nonsingular spacetimes are in- 
compressible Fermi liquids. As a result there is an exact dictionary to the 
quantum Hall system j^Sl IM] along with a Woo algebra of area preserving dif- 
feomorphisms. The grayscale distribution in our analysis describes regions of 
the phase space plane where the lowest Landau level has filling fraction less 
than one. It would be interesting to develop an understanding of this from the 
statistical mechanical point of view, especially with respect to possible relations 
to the fractional quantum Hall effect. 

We have proposed a picture in which the thermodynamic character of very heavy 
states in gravity arises because the effective description in geometry integrates over 
quantum details. As is well known, inertial observers in empty de Sitter space and 
Rindler observers in empty fiat space also see horizons with an entropy and a ther- 
mal character. How might this arise? From our experience in field theory, we are 
used to the decoupling theorem that says that ultraviolet physics decouples from the 
infrared. As a result, integrating out the ultraviolet physics in a field theory vacuum 
state can lead to an infrared effective field theory in a pure vacuum state. However, 
the decoupling theorem fails in theories with gravity. It is tempting to speculate that, 
in this context, integrating ultraviolet details can lead to a low-energy mixed state 
because of entanglement between the infrared and ultraviolet parts of the wavefunc- 
tion of spacetime. Perhaps this entanglement and the resulting mixed nature of the 
low-energy description are more readily visible in some kinematic frames in which an 
observer effectively has large energies with respect to the natural Hamiltonian. 
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A Comments on correlation functions 

A.l Semiclassical distribution functions — another derivation 

Given an arbitrary A^-particle quantum state, when we integrate over — 1 of the 
particles, we obtain a single-particle state described in terms of a density matrix. If 
we consider this procedure among typical states in the ensembles discussed in this 
paper, the average density matrix that we would expect to obtain is 



For example, in the ensemble in which there is no constraint in the number of columns 
of the Young diagram, the collection {c„} is thermal, and precisely given by the Fermi- 
Dirac distribution 



for a suitable chemical potential /i. This is a meaningful statement since we do expect 
the probability of a single fermionic excitation to be n to equal the above distribution. 

As has been argued along Section 4, the average density matrix p^, is related to 
the Husimi phase space distribution as 




(215) 



n>0 



C(„) 



1 + g"-/^ 



(216) 




(217) 



n>0 



where as before z = {p + iq)/\/2h. In the semiclassical limit, where the set of {c„} 
varies very slowly as a function of n, our considerations suggest that 





To show this in more detail, assume that 




(219) 



s 
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with 6„ ^ 1. Then 



Z\ _U|2 



1 

27r^ 



27rn ^ ' 



c(|z|^). (220) 



s 

1 



2'Kh 



Thus, at least in the semiclassical hmit, the phase space distribution function u{p, q) 
should equal the Fermi-Dirac distribution, an statement that we also derived from 
the slope in the limit curve of the limit Young diagram through p56p . 

A. 2 Single trace correlation functions 

We are interested in computing the variance of a correlation function in an ensemble 
of partitions: 

a' = J2^{X){X\F\X){X\F\X) - (^Y.w{X){X\F\X)^ . (221) 

In this expression A corresponds to a partition and we will write n G A to indicate 
that harmonic oscillator level n is occupied when the fermions are in the state cor- 
responding the partition A. As in Section 5, let F be a rotationally invariant sum of 
one-fermion operators. Then 

{X\F\X) = Y,{n\F\n). (222) 

The probability that n G A is 

v{n) = w{X) (223) 

A,neA 

and since we assume that we work in a grand canonical ensemble where the probability 
that n and m are occupied is simply v{n)v{m), we can rewrite cr^ as 



(j^ = w{X){n\F\n){m\F\m) 

A nSA mSA 

— w{X)w{fi){n\F\n){m\F\'m) 

A fJ, neA mefj, 

= v{n){{n\F\n)Y + v{rn)v{n){n\F\n){m\F\m) 
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— f (n)^((?2|F|ri))^ — v{ra)v{n){n\F\n){ni\F\ra) 
= Y^{v{n) - v{nf){{n\F\n)f 

n 

= "^{vin) -v{nf){n\F^\n). (224) 

n 

In the last line we used the rotational invariance of F. This is simply the expecta- 
tion value of in the density matrix Ylni'^i''^) ~ v{n)'^)\n){n\. As we described in 
Section A.l, in the semiclassical limit this becomes an integral over the p, q plane of 
q) ~ "^{Pi times g)^. Notice that for fermionic statistics the variance in 
the occupation number = {ri^) — {n)"^ = (n) — (n)^ since n only takes the values 
and 1. This implies that u — can be interpreted as the variance in the fermion 
occupation number, and the expression here is in fact equivalent to the one given in 
()192j) . Thus, the variance in the expectation value of F is the integral of against the 
variance in the fermion occupation number. In Section 5 we derived a closely related 
expression relating the variance in F to the variance in the distribution function on 
phase space. Following the same reasoning as in that section, one can argue that the 
deviations from universality for delocalized operators goes as . This meets our 
expectations since in the ensembles of interest to us, the entropy scales like N . 
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